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ABSTRACT 


This  thesis  describes  a  circuit  model  for  the  inductive  strip 
in  inhomogeneous  f inline  with  the  following  geometry;  fin  and  strip 
centered  in  the  shield,  dielectric  material  with  e^=2.22, 
b/a  =4/9,  0.5  ^  W/b  ^1.0,  T/a  ^  0.01  and  0.0  ^  d/a  ^0.1  .  The 
model  is  shown  to  produce  results  that  agree  with  data  computed 
using  the  spectral  domain  method.  The  model  has  been  generated 
using  WR(90)  waveguide  operating  in  the  TEjo  mode. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Finline  is  a  transmission  structure  for  electromagnetic 
waves  that  was  first  discussed  by  Meier  in  1974  [Ref.  1]  . 
Finline  consists  of  one  or  more  thin  metal  fins  printed  on  a 
dielectric  substrate  mounted  in  the  E-plane  of  a  rectangular 
waveguide.  Figure  1  depicts  the  particular  variation  under 
consideration  here.  There  are  several  advantages  to  the 
finline  structure.  The  fin  manufacture  is  simplified  by  the 
presence  of  the  dielectric  substrate  which  allows  the  use  of 
well  developed  etching  technologies  for  their  manufacture. 
Finline  has  less  stringent  tolerance  requirements  than 


Inductive  Strip  in  Finline 


^  Dielectric 
_  Metal  Fin  and  Strip 


Figure  1.  Arrangement  of  inductive  strip  in  finline. 


microstrip  and  the  use  of  a  rectangular  waveguide  shield 
simplifies  the  interfaces  with  other  rectangular  waveguides. 
Like  conventional  waveguides,  f inline  can  operate  in  single 
mode  with  low  attenuation.  [Ref.  1] 

Computer  assisted  design  (CAD)  is  essential  to  the 
development  of  complex  devices  in  all  fields  of  engineering. 
This  is  especially  true  in  microwave  and  millimeter  wave 
design.  Numerical  solutions  to  many  electromagnetic  problems 
result  in  time  consuming  computer  implementations .  To  be 
useful  in  design,  a  workstation  should  be  aJble  to  quickly 
predict  the  response  of  a  known  structure.  The  use  of 
equivalent  circuit  models  for  common  structures  that  respond 
in  the  same  way  as  the  general  solution  for  a  limited,  but 
useful,  range  of  parameters  has  been  found  to  be  a  useful  CAD 
technique . 

B.  USE  OF  PREVIOUS  WORK 

Initial  work  by  Knorr  and  Shayda  [Ref.  2]  solved  the 
electromagnetic  fields  in  an  arbitrary  section  of  finline 
using  the  spectral  domain  method.  The  program  IMPED 
implemented  the  solution.  Knorr  and  Deal  solved  for  the 
scattering  parameters  of  an  arbitrarily  located  inductive 
strip  in  finline  and  implemented  the  solution  in  a  program 
called  STRIP  [Ref.  3], [Ref.  4].  Morua  developed  a  circuit 
model  for  homogeneous  finline,  using  input  from  both  IMPED  and 
STRIP  [Ref.  5].  Grohsmeyer  developed  a  model  for 
inhomogeneous  finline  using  the  IMPED  prograun.  This  work  uses 
all  of  the  previous  work  as  a  basis  to  begin  and  relies  on  the 
work  of  Grohsmeyer  for  the  finline  model  [Ref.  6] . 

C .  OBJECTIVE 

The  purpose  of  this  study  is  to  create  a  model  that  allows 
the  representation  of  an  inductive  strip  in  inhomogeneous 
finline.  The  model  must  be  accurate  and  suitable  for 
implementation  in  a  CAD  environment. 
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The  desired  range  of  validity  is  as  follows: 


•  £,  =2.22 

•  W/b  -  ratios  from  0.5  to  1.0 

•  b/a  -  ratio  of  4/9 

•  T/a  -  ratios  greater  than  0.01 

•  d/a  -  ratios  from  0.0  to  0.1 

D .  CONVENTIONS 

Previous  work,  referenced  above,  specified  the  siobscript 
'  eq'  to  denote  an  ecpaivalent  dimension  (length  or  dielectric 
constant) .  This  model  requires  several  equivalent  dimensions 
to  be  specified.  Therefore,  the  following  subscript 
convention  will  be  used. 

•  '  '  -  actual  dimension  (no  subscript) 

•  ' f '  -  equivalent  dimension  in  finline 

•  'd'  -  equivalent  dimension  in  dielectric  loaded  half 

waveguide 

•  'a'  -  dimension  of  the  air  filled  half  waveguide. 


Figure  2  depicts  three  views  of  the  strip  in  finline  and 
labels  the  parameters  used  to  describe  the  geometry.  The 
dotted  lines  are  grid  lines  to  assist  in  specifying  the 
dimensions . 
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•a  -  width  of  the  shield 

•  b  -  height  of  the  shield 

•  T  -  length  of  inductive  strip 

•  d  -  thickness  of  dielectric 
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II.  MODEL  DEVELOPMENT 


A.  CONCEPT 

The  model  for  an  inductive  strip  was  derived  from  the 
homogeneous  inductive  strip  model  developed  by  Morua  [Ref.  5]  . 
The  model  retains  a  similar  structure  but  the  calculation  of 
element  values  has  been  modified  to  account  for  the  inclusion 
of  dielectric. 

The  model  takes  the  physical  structure  of  the  finline 
inductive  strip  and  divides  it  into  two  half-waveguides  that 
are  treated  as  separate  elements,  as  shown  in  Figure  3.  One 


Two  Waveguides 

b 

a/2 

a/2 

Figure  3.  Looking  into  each  of  the  half  waveguides. 


waveguide  is  air  filled  with  transverse  dimension  a/2  and  b. 
The  second  waveguide  is  dielectric  loaded  with  the  sautie 
dimensions  as  the  air  filled  guide.  At  the  mouth  of  each  of 
these  half-waveguides,  there  is  a  parallel  inductor/capacitor 
circuit  that  accounts  for  stored  electric  and  magnetic  energy. 
The  circuits  are  then  connected  through  a  perfect  impedance 
transformer  to  the  f inline.  Figure  4  shows  the  entire 
structure . 


Figure  4.  This  figure  shows  the  circuit  at  one  end  of  the  two 
half-waveguides.  The  other  ends  of  the  waveguides  are 
connected  to  the  mirror  image  of  the  eJDOve  circuit . 


B.  DERIVATION  OF  THE  MODEL 
1 .  General  Comments 

A  model  for  the  loaded  half-waveguide  that  could  be 
represented  in  CAD  software  was  found.  This  required  that  an 
equivalent  below-cutoff  homogeneous  waveguide  be  found.  A 
model  form  that  replicated  the  scattering  properties  of  the 
inductive  strip  as  predicted  by  STRIP  was  then  developed. 

2  .  Loaded  Waveguide  Model 

The  loaded,  below-cutoff  waveguide  was  modeled  by 
matching  the  propagation  constant  and  voltage-power  impedance 
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with  that  of  equivalent  homogeneous  waveguide.  The  first  step 
of  this  process  was  to  write  the  solution  for  the 
electromagnetic  wave  propagation  for  a  below-cutoff  dielectric 
loaded  waveguide  as  shown  in  Figure  3.  In  the  following 
discussion,  it  is  understood  that,  P,  will  be  a  positive 
imaginary  number  due  to  the  below-cutoff  operation  of  the 
guide  and  will  result  in  a  negative  real  propagation  constant. 
The  full  derivation  of  this  solution  is  shown  in  Appendix  A 
and  parallels  the  discussion  on  propagating  loaded  waveguide 
in  Pozar,  [Ref  7;pp.  151-153] .  This  analysis  results  in  the 
following  transcendental  equation 

0=y-^o"P^  tan(d^ej.ico -P^  ) +^ej.ico -P^  tan (  (a-d) ^ico^-P^  )  . 


Since  this  equation  must  be  true  for  all  frecjuencies,  it  is 
possible  to  find  a  set  of.  P's,  for  various  frequencies  and 
geometries.  Several  MATliAB  functions  were  written  to  compute 
P,  in  a  half  section  of  WR(90)  waveguide  for  frequencies 
between  6  GHz  and  the  cutoff  frequency,  and  for  2d/a 
values  from  0  to  0.25  .  Figure  5  shows  the  cutoff  frecjpaency 
as  a  function  of  d/a  for  WR(90)  waveguide.  Figure  6  depicts 
P  versus  frequency  for  various  dielectric  thicknesses .  These 
curves  describe  the  propagation  characteristics  of  a 
particular  geometry  and  are  the  characteristics  that  need  to 
be  matched  by  the  equivalent  waveguide.  The  equation  for  p  in 
a  homogeneous  waveguide  can  be  written  as 


P  = 


(2) 


Using  another  set  of  MATLAB  routines  values  of  a^  and  that 
cause  a  waveguide  described  by  equation  (2)  to  match  the 
curves  in  Figure  6  can  be  found  for  each  desired  geometry. 
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Figure  5.  Cutoff  frequency  vs.  d/a. 


e  7  6  9  10  11  12  13  14 


Frequency  in  GHz 

Figure  6 .  0  vs .  frequency  for  several  values 

of  d/a. 
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That  data  can  then  be  reduced  to  the  following  two  equations 
using  a  least  mean  square  error  curve  fitting  routine; 


6^^=1.0+0.0192  (2-)  -0.3391(2  — )  2+10. 6174  (2  —  )  3  + 

a  a  a 

-3.6647  (2  —  )  <+0.01356  (2-)^ 
a  a 


and 


-^=1.0+0.001314(2-)  -0.03026  (2  — )  ^+0 . 3424  (2  — )  2  + 
a/2  a  a  a  (4) 

-2.1617  (2-)<+1.7461(2  — 
a  a 


where 

Oi  — iO.l. 
a 


The  shapes  of  the  curves  described  by  the  above  equations  are 
shown  in  Figures  7  and  8 .  Application  of  these  equivalent 
dimensions  to  equation  (2)  results  in  p's  that  differ  from 
those  found  for  the  loaded  waveguide  by  less  than  0.2%. 
MATLAB  functions  used  in  the  production  of  the  equivalent 
homogenous  half -waveguide  in  this  section  can  be  found 
Appendix  B . 

The  voltage-power  impedance  of  the  loaded  waveguide 
can  be  calculated  analytically  since  the  electric  and  magnetic 
fields  can  be  found  to  within  a  single  constant  value.  The 
voltage-power  impedance  is  given  by 


(5) 
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In  this  case,  since  the  waveguides  are  operating  below-cutoff 
the  power  will  be  imaginary  and  will  result  in  an  imaginary 
impedance  [Ref.  8] .  The  power  can  be  found  to  be 


.a2a>P^cos2  tan(k^(a-d) 
2kl 


^  (a-d) 

2JcaCos  {k^  ia-d) ) 


A^w^\isin{kad)cos{k^)  ^A^dcjpn) 


2icJ 


2icJ 


(6) 


Equation  (6)  can  be  written  more  compactly  as 

P=b  Cp  (7) 

where  the  imaginary  constant  Cp  represents  the  terms  in 
parentheses  above.  The  voltage  in  equation  (5)  is  given  by 

b 

-fE„dl^-Ej>  (8) 

0 


where  E„  is  the  maximum  electric  field.  The  voltage-power 
impedance  for  the  loaded  half-wave  guide  is 


(9) 


The  voltage-power  impedance  for  the  equivalent  below-cutoff 
waveguide  is  given  by  [Ref.  6] 


(10) 
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Given  that  and  e,.^  are  already  known,  can  be  found  by 
solving  equation  (10)  for  b^  and  substituting  in  from 

equation  (9)  which  results  in 


^  Cp  480nf 


(11) 


Equations  (3)  ,  (4)  and  (11)  define  the  loaded  half-waveguide 
model .  The  variation  in  b^  with  dielectric  thickness  is  shown 
in  Figure  9 . 


Variation  of  b^^with  d/a 


Figure  9.  Variation  in  the  equivalent  height  of  the 
dielectric  loaded  half  waveguide. 
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3 .  Integrated  Model 

Using  the  circuit  model  developed  by  Morua  as  a 
starting  point,  a  MATLAB  function  that  replicated  the 
performance  of  the  CAD  workstation  was  written.  This  program 
was  validated  using  arbitrary  values  of  the  inductance, 
capacitance,  and  turns  ratio  on  the  CAD  workstation.  The  CAD 
program  and  the  MATLAB  functions  used  are  located  in  Appendix 
C.  A  function  minimization  routine  was  then  written  which 
found  values  of  the  various  elements  that  resulted  in  the 
model  output  matching  the  spectral  domain  data  produced  by 
STRIP.  During  the  course  of  much  trial  and  error,  it  was 
found  that  the  expression  used  by  Morua  to  calculate  the  turns 
ratio  in  the  homogeneous  case  could  not  be  extended  to  the 
inhomogeneous  case.  The  model  produced  good  results  when  the 
conductance  on  the  strip  side  of  the  transformer  was  specified 
as  a  function  of  frequency.  It  was  found  that  the  conductance 
could  be  defined  in  the  following  way 

1^0=0 . 001  ( to )  (12) 

where 


Figure  10  depicts  the  way  conductance  varies  with  frequency 
for  a  typical  geometry.  Once  the  value  of  the  conductance 
needed  to  match  the  spectral  domain  data  was  found,  the  turns 
ratio  of  the  impedance  transformer  was  calculated  using  the 
relation 


(14) 
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where  N  is  the  turns  ratio,  Y,  is  the  finline  conductance,  Yg 
is  the  conductance  on  the  strip  side  of  the  transformer.  The 
behavior  of  the  turns  ratio  with  frequency  is  shown  in  Figure 
11.  This  model  was  found  to  be  able  to  replicate  the 
performance  of  the  inductive  strip  for  all  geometries  of 
interest  with  less  than  1%  error.  Some  results  are  shown  in 
Figures  12  and  13  (in  the  figures  spectral  domain  data  is 
referred  to  as  SPEC_DOM  and  the  model  is  referred  to  as 
STRIP)  . 


x10“^  Conductance  Variation  with  Frequency 


Frequency  (GHz)  W/b-0.5,  T-50,  and  d—30 


Figure  10.  Variation  of  the  conductance  with  frequency. 


Turns  Ratio 


Variation  in  Turns  Ratio  with  Frequency 


Figur*  11.  Variation  in  Turns  Ratio  with  Frecjuency. 
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12 .  Response  of  Sn  for  Model  and  Spectral  domain  data 
for  M/b*0.5,  T»10  mils,  d»0  mils,  L»14.67nH,  C-0.0043pF  and 
Yo-0 . 001  (1 . 315+3 . 63f„-3 . 43f„*+0 . 98f.’)  . 
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rigur*  13.  The  response  of  for  Model  and  Spectral  domain 
data,  forW/b«0.5,  T«50  mils,  d«30  mils,  L«12.12nH,  C*.0018pF, 
and  Yo*0. 001(1. 69+2. 06f.-1.52f^+0.3792f„’)  . 


PHASE 


IZI.  DATA  COLLBCTZON 


A.  STRIP  PROGRAM  DBSCRZPTIOH 

The  input  data  was  generated  using  the  STRIP  program. 
STRIP  is  a  Fortran  prograun  that  computes  the  scattering 
pareuaeters  of  the  inductive  strip  in  f  inline  using  the 
spectral  domain  method.  The  STRIP  prograun  and  the  development 
of  the  theory  behind  it  are  discussed  in  detail  in  References 
3  and  4 .  Accuracy  of  1-2%  with  experimental  data  on 
homogeneous  finline  were  reported.  [Ref.  3], [Ref.  4] 

B.  NODiriCATIOMS  TO  THB  ORIGINAL  PROGRAM 

During  the  course  of  the  data  collection  both  a  SUN 
workstation  and  a  VAX  11/785  were  available  to  compute  the 
scattering  data.  The  two  computers  were  found  to  produce 
results  that  differed  by  as  much  as  10%.  These  errors  were 
presumed  to  result  from  differences  in  the  way  the  two 
machines  handled  small  numbers.  This  indicated  a  possible 
problem  in  the  way  the  program  dealt  with  small  numbers.  To 
correct  the  problem,  the  program  was  cheuiged  to  re-specify 
some  variables  that  were  previously  real  as  double  precision. 
This  modification  elevated  the  problem  and  resulted  in 
agreement  between  the  solutions  generated  by  the  VAX  and  the 
SUN.  Unfortunately,  other  intermittent  problems  began  to 
appear  in  the  output  data.  Further  extensive  trovible  shooting 
revealed  that  some  double  precision  values  were  being  returned 
from  subroutines  into  real  nvimbers .  The  extra  bits  were 
written  into  the  adjacent  memory  spaces,  corrupting  them. 
This  problem  was  corrected  by  re-compiling  the  progreun  to 
treat  all  real  numbers  as  double  precision.  This  final 
alteration  corrected  all  the  apparent  problems  in  the 
operation  of  the  program. 
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C .  DATA  GATHERING 

The  data  on  which  the  model  was  based,  was  calculated 
using  parameters  associated  with  WR(90)  waveguide.  WR(90) 
waveguide  was  selected  for  two  reasons.  First,  it  is  more 
convenient  to  conduct  measurements  in  X-band.  Second,  as 
shown  by  Morua  in  the  homogeneous  case,  selection  of  b/a=4/9 
may  permit  the  model  to  be  extended  to  b/a's,  from  0.4  to  0.5. 
The  following  list  specifies  the  values  of  the  parameters  used 
in  determining  the  data  points .  All  measurements  in  mils 
(1  mil  =  0.001  inches): 

•  W/b  =  {  0.5,  0.7,  0.9,  1.0  } 

•  T  =  {  10,  20,  30,  40,  50,  60,  70,  80,  90,  100,  200,  300, 

400,  500,  600,  700,  800,  1000  } 

•  d  *  {  0.0,  10,  30,  50,  70,  90} 

The  strip  length  interval  was  increased  for  T>100  mils  because 
the  scattering  par2uneters  were  found  to  vary  more  slowly  at 
large  values  of  T. 

The  performance  of  the  STRIP  progrzun  is  determined  in 
large  part  by  the  order  of  the  matrix  of  inner  product  terms. 
Morua  showed  that  most  of  the  convergence  for  typical 
parameter  values  had  occurred  by  matrix  of  order  10. 
Therefore,  a  matrix  of  order  10  was  used  as  a  starting  point 
to  calculate  the  data .  Some  data  sets  were  re-computed  with 
matrix  orders  of  12  to  attempt  to  correct  errors  discussed 
below. 

The  STRIP  prograun  was  developed  based  on  the  assumption 
that  the  f inline  was  operated  in  the  TE^q  mode  and  the  two 
half -waveguide  sections  were  operated  below  cutoff.  As  the 
thickness  of  the  dielectric  substrate  was  increased  the  second 
propagating  mode  began  to  have  an  effect  for  the  higher 
frequencies  in  the  range  of  interest.  In  this  situation,  the 
program  returns  an  increasing  S^  magnitude  normally  1.0. 
Values  of  the  magnitude  of  Su  were  assumed  to  be  continually 
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decreasing  with  increasing  frequency.  This  assumption 
complies  with  intuition  of  the  physical  system  that  as  the 
non-propagating  half-waveguide  approaches  the  conditions  for 
propagation,  more  energy  will  be  transmitted  and  less 
reflected.  Therefore,  data  points  that  indicated  an 
increasing  magnitude  of  with  frequency  were  discarded.  At 
apparently  random  geometries,  the  STRIP  program  returned 
values  for  the  magnitude  of  Sn  that  did  not  correspond  with 
the  similar,  but  different,  geometries  or  the  expected 
behavior  of  the  structure.  Identification  of  this  data  is 
discussed  below.  The  block  of  erroneous  data  would  typically 
begin  at  a  strip  length,  T,  of  approximately  500  mils  and 
continue  through  1000  mils  for  a  given  dielectric  thickness 
and  gap  width.  This  type  of  error  was  found  to  effect  15%  of 
the  data.  By  not  including  data  with  T  greater  than  500  mils 
the  amount  of  affected  data  was  reduced  to  approximately  5%. 
For  this  reason,  data  with  a  strip  length  greater  than 
500  mils  was  discarded.  An  exeimple  of  this  behavior  is  shown 
in  Teible  1  (on  the  following  page)  .  The  bad  data  is  located 
at  d*70  mils,  T*=600  mils  and  f=12  GHz.  The  numbers  in  the 
table  represent  the  magnitude  of  S^.  The  physical  arrangement 
of  the  system  is  such  that  with  increasing  strip  length,  T, 
the  reflection  coefficient  should  continually  increase.  The 
adjacent  data,  at  d=50,  does  exhibit  the  expected  behavior. 
The  data  point  described  here  and  others  that  could  be 
identified  as  bad  in  this  manner  were  discarded. 
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Table  1.  Computed  Magnitude  of  near  d=70  and  Ts€00 
at  12  GBz,  Illustrating  Anomalous  Results 


f=12GHz 

d=50  mils 

d=70  mils 

T=300  mils 

0.8790 

0.7877 

T=400  mils 

0.9643 

0.8216 

T=500  mils 

0.9788 

0.8242 

T=600  mils 

0.9872 

0.6259 

21 


III.  SLSMENT  REDUCTION 

A.  CONCEPT 

It  has  been  shown  that  element  values  could  be  found  which 
produced  excellent  agreement  with  the  STRIP  data  for  specific 
geometries .  The  development  of  a  useful  tool  for  the  CAD 
environment  requires  that  a  method  of  generating  element 
values  as  a  function  of  the  geometry  be  found.  In  the 
following  discussions,  the  coefficients  of  the  strip 
conductance  are  referred  to  as  elements . 

B .  REDUCTION  TECHNIQUE 

Polynomials  were  selected  as  generating  equations. 
Polynomials  were  selected  because  of  their  simplicity  and 
because  the  element  data  appeared  to  be  slowly  varying.  The 
function  minimization  routines  resulted  in  solutions  that  were 
exact  beyond  the  level  needed  for  a  successful  model.  Errors 
of  less  than  1%  with  the  STRIP  data  were  typical  and  the 
expected  error  for  the  STRIP  data  was  1~2%.  Therefore,  it  was 
assumed  that  an  acceptable  solution  could  be  found  and  that  a 
least  mean  square  curve  fit  would  be  an  adequate  method  of 
reducing  the  data.  Progreuns  were  written  that  conducted  a 
least  mean  square  fit  over  each  of  the  three  parameters  of 
interest . 

The  reducing  function  first  took  an  entire  set  of  data  for 
a  given  W/b  and  then  subdivided  that  into  sets  of  different 
d/a  values.  For  each  set  of  like  d/a's  a  second  order  least 
mean  square  reduction  was  done  using  the  T/a  ratios  as  the 
input  parameter.  The  resulting  coefficients  were  placed  in 
column  form.  This  procedure  was  repeated  for  each  different 
d/ a  set.  The  coefficients  were  then  arranged  into  data  sets 
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treating  each  location  in  the  colxamn  as  a  data  point .  Another 
least  mean  square  fit  was  done  using  the  d/a  ratios  as  the 
input  pareuneter.  Another  set  of  coefficients  were  generated 
from  this  reduction.  These  coefficients  were  again  placed  in 
a  known  column  order.  This  procedure  was  repeated  for  each 
W/b  set.  The  coefficients  were  again  placed  in  a  known  column 
order  and  the  data  reduced  to  coefficients  using  a  least  mean 
square  fit  with  W/b  as  the  input  parameter.  This  procedure 
resulted  in  162  coefficients.  The  element  values  could  then 
be  regenerated  by  reversing  this  process .  The  progreuns  that 
implemented  the  above  algorithm  were  tested  by  picking  an 
arbitrary  set  of  coefficients,  generating  element  values  then 
reducing  the  element  values  to  coefficients.  The  coefficient 
values  were  recovered  with  differences  on  the  order  of  10"®. 
This  test,  of  course,  did  not  give  any  indication  of  the 
performance  of  the  algorithm  with  unsmoothed  element  values . 
The  function  listings  for  the  above  routines  are  included  in 
Appendix  C. 

Application  of  the  procedure  described  above  to  the 
element  values  was  unsatisfactory.  Testing  the  model  with 
elements  generated  from  the  curve  fitting  routine  resulted  in 
large  errors .  Typical  error  levels  were  found  to  be  in  the 
neighborhood  of  80%.  Apparently,  small  errors  generated  in 
the  initial  recursions  were  magnified  by  using  curve  fitting 
on  the  coefficient  in  the  subsequent  curve  fits.  After  the 
initial  attempt  was  found  to  be  inadequate,  several  additional 
procedures  were  used  to  attempt  to  develop  suitable 
coefficients . 

The  element  values  were  examined  and  values  that  were  far 
from  the  apparent  trend  of  the  data  were  eliminated  to  create 
a  data  set  that  could  be  used  to  generate  usable  coefficients. 
This  procedure  did  improve  the  performance  but  the  results 
were  again  found  to  be  unsatisfactory.  Several  possible 
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variations  in  the  'correct'  trend  were  explored  with  similar 
unsatisfactory  results. 

The  most  promising  of  the  coefficient  sets  generated  zdsove 
was  used  to  generate  all  of  the  values  with  the  exception  of 
the  value  of  the  inductor.  The  inductor  was  selected  because 
it  was  the  single  most  significant  element  to  the  model's  fit. 
Another  program  was  written  that  found  the  value  of  the 
inductor  needed,  given  that  the  other  elements  were  all 
generated  by  coefficients.  The  elements  were  again  reduced 
using  the  initial  procedure  described.  A  comparison  between 
desired  element  values  and  those  generated  by  the  coefficients 
was  made  and  showed  that  the  elements  other  than  the  inductor 
were  a  perfect  fit,  as  expected.  The  error  in  the  value  of 
the  inductor  was  found  to  be  in  the  5-7%  range,  with  some  much 
larger  values .  This  level  of  error  in  the  inductor  was  found 
to  result  in  20-30%  error  when  compared  to  the  STRIP  data. 
This  procedure  was  applied  successively  to  each  element  in  the 
hope  that  the  coefficient  values  would  converge  to  an 
acceptable  solution.  Convergence  did  not  occur.  The  output 
of  the  data  was  continually  found  to  generate  coefficients 
that  produced  error  levels  in  the  20-30%  range.  Appendix  E 
contains  a  portion  of  the  MATXiAB  output,  with  typical  output 
values . 

Two  final  attempts  at  function  minimization  were  made. 
Both  attempts  iterated  through  the  entire  data  set  to  allow  a 
best  guess  at  the  correct  coefficients.  The  first  attempt 
used  the  coefficients  found  using  the  curve  fitting  procedure 
to  generate  all  of  the  element  values  except  the  value  of  the 
inductor.  The  27  coefficients  necessary  to  specify  the 
inductor  were  passed  to  the  function  minimization  routine . 
This  program  tested  every  fifth  data  point  of  the  entire  STRIP 
data  set.  The  error  at  each  point  was  defined  in  the  same  way 
as  described  in  the  previous  chapter .  The  total  error  was 
defined  as  the  summation  of  the  error  at  all  the  test  points . 
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The  resulting  coefficients  produced  no  significant  change  in 
the  performance  of  the  model.  A  modification  to  the  program 
was  made  that  allowed  the  function  minimization  routine  to 
vary  all  162  of  the  coefficients.  This  modification  did  not 
result  in  improved  performance.  The  failure  of  these  last  two 
attempts  was  expected  for  two  reasons.  First,  the  function 
minimization  routine  was  designed  to  operate  with  less  than 
five  parameters.  Second,  the  nature  of  the  model  is  such 
that  when  drastically  incorrect  element  values  are  applied  to 
the  model,  the  value  of  the  magnitude  of,  S^,  goes  to  either 
0  or  1  and  the  angle  to  ±180*,  as  a  result,  no  indication  of 
increased  or  decreased  error  is  given  to  the  minimization 
routine . 

Coefficients  were  not  found  which  allowed  the  element 
values  to  be  represented  as  series  of  equations .  Possible 
alternative  approaches  are  discussed  in  the  following  chapter. 
The  reason  that  least  mean  square  curve  fitting  was  not 
applicable  is  not  completely  clear,  but  it  is  most  likely  a 
combination  of  the  following  factors. 

•  The  output  of  the  STRIP  program  changed  after  the 
variables  were  re-specified  as  double  precision. 
Indicating  a  possible  problem  in  the  program,  although  the 
program  has  been  confirmed  experimentally  for  several 
geometries . 

•  The  STRIP  data  may  contain  unlocated  irregularities  that 
prevent  the  production  of  smooth  element  values. 

•  Function  minimization  accentuates  inconsistencies  in  the 
data  by  finding  the  optimum  fit  point  by  point. 

•  Using  an  element  generation  method  that  results  from  three 
iterations  of  curve  fitting  is  susceptible  to  significant 
error  magnification  during  the  calculation  of  the 
elements . 
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V.  SUMMARY 


A.  CONCLUSIONS 

The  equivalent  circuit  model  presented  has  been  shown  to 
accurately  reproduce  the  response  predicted  by  the  STRIP 
prograun  for  individual  geometries.  Within  the  model,  there 
are  three  effects  that  act  together  to  replicate  the  behavior 
of  the  actual  inductive  strip.  They  are  the  decay  of  the 
field  along  the  length  of  the  inductive  strip,  the  scattering 
from  the  edge  of  the  strip,  and  the  coupling  of  the  field 
between  the  f inline  and  the  half -waveguide  sections.  The 
length  of  the  inductive  strip  was  accounted  for  by  the  two 
below-cutoff  waveguides,  which  attenuated  the  energy  flow 
along  their  lengths .  The  scattering  caused  by  the  edge  of  the 
inductive  strip  was  modeled  with  the  fixed  parallel  inductor 
and  capacitor  circuit.  The  coupling  of  the  field  from  the 
finline  to  the  below-cutoff  waveguides  was  modeled  with  a 
frequency  dependent  conductance  defined  at  the  beginning  of 
the  strip.  This  frequency  dependant  conductance  accounted  for 
the  transformation  of  the  finline  conductance  to  the 
conductance  of  the  below-cutoff  waveguides .  Viewing  the  model 
in  this  way,  gives  valuable  insight  into  the  behavior  of  the 
actual  inductive  strip  in  finline. 

B .  RBCCBOOCNDATIONS 

An  algorithm  which  accurately  reproduced  the  element 
values  was  not  found.  To  be  useful  in  the  CAD  environment,  a 
method  of  generating  the  element  values  from  the  geometry  is 
essential.  One  possibility  would  be  to  smooth  the  STRIP  data 
using  some  curve  fitting  technique  prior  to  applying  the  model 
to  the  data.  Smoothing  the  data  first  would  reduce  the 
tendency  of  the  function  minimization  routines  to  produce 
element  values  that  are  difficult  to  describe  analytically. 
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Element  values  that  can  be  easily  and  accurately  reduced  to 
equations  are  essential  because  of  the  three  levels  of 
generation  required. 

The  model  should  be  extended  to  other  geometries .  It  has 
been  tested  at  values  of  W/b  as  low  as  0.1  and  has  continued 
to  perform  well  for  specific  geometries .  The  model  may  be 
scalable  to  b/a's  of  0.4  to  0.5  as  was  the  case  for 
homogeneous  finline,  but  this  will  require  further 
investigation . 

The  model  for  the  loaded  waveguide  could  be  tested 
experimentally  and  implemented  directly  in  CAD  software. 
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APPENDIX  A.  DERIVATION  OF  LOADED  HALF  WAVEGUIDE  MODEL 


This  appendix  presents  a  derivation  of  the  attenuation 
through  a  partially  loaded  below  cutoff  waveguide  and  the 
computation  of  the  equivalent  height  by  matching  the  voltage- 
power  impedance.  The  variable  'a'  will  be  used  here  to  denote 
the  width  of  the  shield.  The  presentation  here  follows 
Pozar's  development  for  the  loaded  waveguide.  [Ref  7;pp.  151- 
153] 


As  can  be  seen  in  Figure  14,  the  geometry  is  uniform  in  the  y 
direction  and  since,  n=0,  the  TEio  modes  have  no  y  dependence 
the  wave  equation  for  H,  can  be  written  separately  as 

(1) 
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and 


(2) 


where  k.  and  k^  are  the  wave  numbers  in  the  two  regions  (a  in 
air  and  d  in  the  dielectric) .  These  two  equations  are  related 
by  the  fact  that  both  waves  must  propagate  down  the  waveguide 
together.  Therefore,  the  propagation  constant,  y,  must  be  the 
Seune  on  each  side.  The  model  assumes  no  loss,  therefore,  the 
propagation  constant  reduces  to  j^.  Since,  the  model  also 
assumes  below  cutoff  operation  3  must  be  a  positive  imaginary 
number,  which  will  result  in  a  negative  real  propagation 
constant  and  evanescent  electromagnetic  waves.  Therefore,  3 
can  be  written  as 


Given  the  wave  equations  above  and  a  knowledge  of  the  boundary 
conditions,  the  general  solution  can  be  written  as 

».d  =Acos  (k^)  +flsin(Jc^)  for  x  <  d  (4) 

Hg^=Ccos  {kf^{a-x) ) +Dsin(k^{a-x) )  for  d  <  x  <  a. 


The  field  components  of  a  TE  field  can  be  found  from 


„  _  ju)\x 

"■  jcl 


(5) 


The  electric  field  in  the  y  direction  can  be  found  to  be 

^y<r  (flcos  (kj  x)  -Asin(kj  x) )  (6) 

and 
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{Dcos  (iCg (x-d) )  -Csin ik^id-x) )  . 


(7) 


E=^ 
y  L- 


Boundary  conditions  on  the  tangential  electric  field  imply 
that  B=D=0 .  The  field  components  reduce  to 


Hg^=Acos  (k^  x) 


/4a  =  Ceos  (iCa  (a-x)  ) 


(8) 


Eya=-^^^sin(ka  x) 


Eya= Sin  (kja-x))  (9) 


Sin  ik^) 


/4a=— i:5^Sin(iCa(d-X)  )  .  (10) 


Since,  the  electric  field  tangential  to  the  dielectric 
boundary  must  be  continuous,  the  following  condition  applies 

^,^(d)=E,a(d)  .  (11) 


Similarly,  the  tangential  magnetic  field  must  be  continuous 
across  the  boundary 


(12) 

Eq[uations 

following 

(12)  and  (13)  can  be  solved  and  result 

transcendental  equation 

in  the 

0=iCatan (kjjd)  +k^tan{k^[a-d) )  . 

(13) 

Note  also  that 

cos(jcjd) 
cos  {k^{a-d)  ) 


(14) 


Solving  equation  (3)  for  k,  and  k^  then  substituting  into 
equation  (14)  produces 
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0=^fkf^  tan  (d^e^ico -P^  ) tan  { (a- d)  yfkf^  )  . 


Solutions  for  the  propagation  constant  can  then  be  found 
numerically.  For  a  homogeneous  waveguide  operating  in  the  TE^o 
mode  P  can  be  found  to  be 


(16) 


Equation  (15)  will  produce  a  set  of  P's  for  different 
frequencies  and  geometries .  By  selecting  a^  and  e^d  in 
equation  (18)  correctly,  a  homogeneous  waveguide  can  be  found 
that  approximates  P's  frequency  variation  for  a  specific 
geometry . 

The  voltage-power  impedance,  is  defined  as 


(17) 


Since,  this  is  a  below-cutoff  waveguide,  no  real  power  will 
propagate.  The  imaginary  voltage-power  impedance  can  still  be 
defined  by  using  the  complex  power  [Ref.  8].  The  complex 
power  can  be  written  as 


J»=l 


J  j" J 

oVo  d  , 


dy. 


(18) 


There  is  no  y  variation,  so  the  integration  in  y  can  be  done 
immediately 


P=h 


a 

>0  d 


(19) 
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Equation  (21)  can  be  more  compactly  written  as 

i»=i>Cp  (22) 

where  Cp  is  the  value  of  the  terms  in  parenthesis  in  ecpaation 
(21)  .  Taking  into  account  that  the  powor  is  imaginary  the 
average  power  can  be  written  as 

P=-IM(P)  .  (23) 

2 
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The  voltage  is  given  by 


where  E,  is  the  max  field  given  found  by  finding  the  zero  of 
the  derivative  of  the  y  component  of  the  electric  field  with 
respect  to  x 


dbc 


_ -j-Aojicos  (ic^)  sin(/c,  (x-a) )  ^ 
k^co9{k^{a-d) )  J ‘ 


This  results  in 


O=cos  (ic,  (x-a)  )  —  x=(n+-i)-^+a  n=.  .  . -l,  0,  l .  .  (26) 

2  k^ 


which  gives 


-jAwjiCOS  (iCjd) 

k^cos(k^{a-d) ) 


The  voltage-power  impedance,  can  then  be  written  as 


z  =_X^=jb£:MJ- 

^ov  J  p  p 

^^ave  '-p 


The  voltage-power  impedance  for  the  equivalent  below  cutoff 
waveguide  can  then  be  found  as  [Ref.  8] : 


Z^^=120n 
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Solving  equation  (29)  for  and  substituting  in  the  voltage- 
power  impedance  found  in  equation  (28) ,  defines  the  equivalent 
height  required  to  produce  the  szime  voltage-power  impedances 
for  bot  •  -“guides 

_  PO) 

^  Cp  480iif 
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APPENDIX  B.  MATLAB  FUNCTIONS  FOR  HALF  WAVEGUIDE  MODEL 


This  appendix  presents  an  explanation  and  listing  of  the  MATLAB 
functions  used  to  determine  the  model  for  the  half  waveguide. 
VarieUbles  defined  as  global_' name'  are  global  variables  which  must 
be  specified  in  the  MATLAB  shell.  The  purpose  of  these  functions 
is  to  produce  a  set  of  equivalent  a^'  s  and  s  that  approximate 
the  variation  of  P  for  a  specific  geometry. 

The  cutoff  frequency  was  found  using  FIND_FC.M.  FIND_FC.M  uses 
the  MATLAB  function  minimization  routine  to  find  the  zeros  of  the 
transcendental  equation  contained  in  TRANS_FC.M. 

The  equivalent  a^' s  and  Era's  were  found  using  MODEL_EQ.M.  This 
function  found  the  set  of  p's  over  the  freqpiency  of  interest  for  a 
particular  geometry  by  finding  the  zero's  of  the  transcendental 
function  found  in  Appendix  A  (equation  15)  .  Equation  (15)  was 
implemented  in  TRANS. M.  The  equivalent  dimensions  were  found  by 
matching  the  P's  just  produced  with  P's  calculated  using  test 
dimensions  supplied  by  the  function  minimization  routine.  The 
routine  selected  the  best  choice  of  a^' s  and  Era's  to  fit  the  data. 
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The  following  functions  relate  to  cutoff  frequency  calculations 


function  fc=f ind_fc () 

%  Finds  the  cutoff  frequency  for  the  loaded  half  waveguide. 

%  The  varieUole  da  refers  to  d/a  and  fc  is  the  cutoff  frequency. 
%  The  following  variables  must  be  declared  global .  global  da 

tol=1.0e-5;  %  sets  the  accuracy  of  the  solution 

for  da= . 01 : . 01 : . 25  %  iterates  through  da 

global_da=da;  %  sets  to  global,  pass  to  trans_fc 

freq=13e9;  %  start  of  search 

freq=fzero (' trans_fc' , freq, tol) %  finds  solution  to  trans_fc 
fc=[fc;  [  da  freq/le9]]%  stores  the  solution 
end 


function  y=trans_f c (x) 

%  finds  the  cutoff  frequency. 

beta=0;  %  condition  for  cutoff 

a=. 02286/2;  %  wr90  in  meters 

d=global_da*a;  %  d  dimensions 

er=2 .22; 

ko=2*pi*x/3e8;  %  wave  numbers 

kd=(er*ko''2+beta''2)  ''O.S; 
ka=(ko''2-»-beta''2)  ''0.5; 

%  transcendental  equation  ->  0 
global_yy=ka*tan (kd*d) +kd*tan (ka* (a-d) )  ; 
y=global_yy; 
end 
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The  following  functions  relate  to  the  computation  of  the  equivalent 
dimensions : 

function  y=model_eq (pass _ tol) 

%  global  varieibles  that  must  be  declared:  global_pass_tol, 

%  global_aeq,  global_er,  global_f,  global_y_error  global_yy 
%  g_cost 

global_pass_tol=pass _ ^tol;  %  sets  tolerance  throughout  routine 

tol=pass _ ^tol; 

freq=12;  %  upper  bound  of  matching 

ii=l;  X (ii, :)=[. 01111  1.0701];%  starting  point  for  search 
for  i= . 01 : . 01 : . 25  %  iterates  through  geometries 

n=l; 

hhold=245;  %  beginning  point  for  beta 

global_d=i; 

for  f »6 ; 0 . 5 ; f req  %  iterates  through  frequency 
global_f=f *le9;  %  pass  frequency  to  global 

%  get  beta  at  that  freq 
hhold=f zero ( ' trans' , hhold, le-6) ; 
g_cost (n, 1) =hhold; %  save  the  data  point 
n=n+l; 
end 

%  find  the  aeq  and  er  that  match 
X (ii, : ) =fmins (' cost_eq' , x (ii, : ) , tol) ' ; 

%  the  beta's  held  in  the  cost 
%  matrix  below, 

%  save  data  and  error 

y2(ii,:)=[i  x(ii,:)  global_yy  global_y_error  ] ; 
ii=ii+l ; 

X (ii, : ) =x (ii-1, : ) ; 
end 
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save  y2 

y=y2 

end 


function  y=»trans  (x) 

%  This  function  finds  the  degree  to  which  the  transcendental 
%  equation  goes  to  0  for  a  given  beta  and  dimension. 

beta=x;  %  the  value  being  found 

%  constants 

a*. 02286/2;  d=global_d*a;  er=2.22; 
ko=2*pi*global_f/3e8; %  wave  numbers 
kd=  (er*ko''2+beta‘'2)  ''0.5; 
ka=  (ko''2+beta''2)  ''0 .5; 

%  transcendental  equation 
global_y=ka*tan (kd*d) +kd*tan (ka* (a-d) )  ; 
y»global_y; 
end 


function  y=cost_eq(x) 

global_aeq=x  (1)  ;  %  variedjles 

global_er=x (2)  ; 

tol=global_pasa_tol;  %  set  tolerances 
freq®12;  %  set  upper  limit  on  matching 

hhold=245;  %  starting  value  for  search 

n*!; 

for  f=6:0.5;freq 

global_f*f *le9;  %  passes  frequency  to  trans_eq.m 

%  finds  the  beta  that  solves  for 
%  a  particular  aeq  and  er 
hhold=f zero ( ' trans  eq' , hhold, tol) ; 
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g_cost (n, 2) =hhold; %  saves  data  in  matrix 
n=n+l; 
end 
yl=0; 

for  i=l;n-l 

%  compares  the  beta  vectors  from 
%  load  waveguide  with  homogeneous 
%  waveguide,  error  is  the  square 
%  of  the  sum  of  the  differences . 
yl=  (g_cost  (i,  1)  -g_cost  (i,  2)  )  ''2+yl; 
end 

y_error=yl ; 

y=yl; 

end 


function  y=trans_eq (x) 

%  function  returns  the  beta  for  a  specific  homogeneous  geometry 
%  at  a  specific  frequency 

beta=x;  %  desired  variable 

ko=2*pi*global_f/3e8  %  calculation  of  wave  numbers 
lca=  (global_er*ko''2+beta''2)  ''0.5; 

yy=tan (ka* (global_aeq) ) ;  %  goes  to  zero  to  satisfy  B.C. 

y=yy; 

end 
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APPENDIX  C.  EBSOF  HODEL  AND  KATLAB  MODEL 


This  appendix  describes  the  steps  taken  to  develop  the 
circuit  model.  The  first  step  in  the  process  was  writing  a 
MATLAB  progreun  that  could  replicate  the  performance  of  the  CAD 
software.  The  CAD  software  availed>le,  and  the  target  of  this 
work,  was  the  Touchstone  software  by  EESOF .  Following  this 
discussion  is  a  printout  of  the  structure  described  in  EESOF 
code,  a  typical  EESOF  data  file  and  the  MATLAB  functions 
mentioned. 

The  approach  taken  to  develop  the  MATLAB  code  was  to 
represent  each  of  the  half  waveguides  as  a  7C-equivalent 
circuit,  as  shown  in  Figure  15  (the  dotted  lines  represent 
each  half  waveguide  element) .  The  nature  of  the  structure 
lends  itself  most  easily  to  use  conductance  for  calculations. 
The  value  of  the  individual  elements  of  the  7i-equivalent 
circuit  for  a  waveguide  can  be  found  as 


sinh(Yat) 

V 

tanh(Yat/2) 

(1) 

Y  —  ^ 

sinh(Yjt) 

tanh(Ydt/2) 

where  the  subscript  '  d'  and  '  a.'  refer  to  the  loaded  and 
unloaded  half-waveguides .  The  initial  implementation  of  the 
aJDOve  algorithm  was  done  in  a  function  called  STRUCTURE. M. 
This  function  and  the  seune  model  design  implemented  in 
Touchstone  was  tested  over  a  wide  range  of  parzuneter  values  to 
ensure  that  the  two  programs  produced  the  saone  results.  In 
the  interest  of  speed  STRUCTURE. M  was  divided  into  two 
functions  STRUCTURE_1 . M  and  STRUCTURE_2 . M .  STRUCTURE_2 .M 
computed  values  that  were  dependant  only  on  geometry  and 
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7c- Equivalent  Circuit  Representation 


rigur*  15.  The  rt-equivalent  circuit  for  the  two  half 
waveguides .  The  dotted  lines  separate  each  of  the 
waveguides . 

therefore,  could  be  saved  while  the  values  in  STRUCTURE_2 . M 
depended  on  the  element  values  and  needed  to  be  recalculated 
for  each  iteration. 

The  '  structure'  functions  were  then  called  by  a  series  of 
function  minimization  routines.  Initially  FIRST_SEARCH .M  was 
used  to  get  a  rough  idea  of  the  progression  of  the  element 
values.  This  prograun  first  searched  over  the  coefficients  of 
the,  Yo,  polynomial  to  get  a  solution  reasonably  close  to  the 
values  from  STRIP.  Then,  another  function  minimization  was 
done  using  all  of  the  parameters.  The  purpose  of  this  search 
regime  was  two  fold:  a  search  over  six  parameter  is  very  time 
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consuming  and  the  results  from  one  geometry  to  the  next  can  be 
very  discontinuous.  By  holding  the  value  of  the  inductor  and 
the  capacitor  fixed  for  the  first  search,  the  entire  progrzun 
ran  faster  and  produced  smoother  results.  After  returning  the 
'best'  guess  solution  of  the  parauneters  to  the  calling 
function,  the  pareuneter  values  were  then  stored  and  passed  to 
the  next  iteration  as  the  starting  point  of  the  search.  The 
function  STORE. M  opened  a  file  named  STORAGE. MAT  and  save  the 
results  in  a  column  form.  Figure  16  contains  a  flow  graph 
which  summarizes  the  cdDove  discussion. 

The  performance  of  a  function  minimization  is  determined 
by  the  way  the  error  is  defined.  The  error,  in  this  case,  was 


Figure  16.  The  iteration  between  the  functions  to  produce 
element  values  can  be  seen  above . 
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calculated,  in  the  function  CALC_ERROR.M  using  the  square  of 
the  svun  of  the  differences  of  the  calculated  values  of  the 
scattering  parameters  and  the  values  generated  by  the  STRIP 
progreun.  Because  both  the  phase  and  magnitude  of,  were 
considered  equally  important,  the  square  of  the  difference  of 
the  two  values  was  also  taken  and  added  to  the  error  to  force 
the  errors  to  become  evenly  balanced.  The  following  is  an 
algebraic  description  error 

(|au.l-|s„.|)^  <2> 

(3) 

where,  S^,,  is  a  vector  that  contains  values  of  the  spectral 
domain  data  over  the  frequency  set  and,  Snof  is  a  vector  that 
contains  the  values  calculated  by  the  model .  The  errors  were 
combined  in  the  following  way 

®tocaJ 

After  several  runs  of  FIRST_SEARCH  the  values  of  elements 
were  smooth  enough  to  allow  the  calculation  of  equation  to 
compute  the  element  vales  at  specific  geometries.  The  function 
MOD_ELEMENTS .M  allowed  the  element  values  to  be  hand  tuned  to 
manually  smooth  the  data  prior  to  attempting  to  reduce  the 
data.  The  coefficients  were  calculated  using  FIND_COEF.M  and 
FIND_WB.M.  These  functions  were  placed  into  another  function 
called  MAKE_COEF.M  This  program  iteratively  called  up  the 
saved  values  of  the  elements  and  computed  coefficients  using 
a  least  mean  square  curve  fitting  routine.  For  a  given,  d/a, 
and,  W/b,  the  second  order  polynomial  for  the,  T/a,  variation 
was  found  for  each  of  the  elements.  These  coefficients  were 
then  placed  in  column  form.  The  set  of  data  for  the  next, 
d/a,  geometry  was  reduced  until  all  the,  d/a's,  for  a 
specific,  W/b,  geometry  was  computed.  These  coefficients  were 
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then  treated  as  elements  and  converted  into  polynomials  that 
depended  on,  d/a.  These  coefficients  are  placed  in  colximn 
form  and  the  data  for  each,  W/b,  was  then  reduced.  These 
coefficients  were  again  treated  as  elements  and  polynomials 
were  formed  using  the,  W/b,  dependence.  By  reversing  the 
process  in  the  functions  GETJWB.M  and  GET_COEF.M,  the  values 
of  the  elements  were  regenerated.  This  program  was  tested 
extensively  by  starting  with  an  arbitrary  set  of  coefficients 
and  generating  elements  then  re-created  the  coefficients  and 
verifying  that  the  original  coefficients  were  regenerated. 

After  producing  a  set  of  coefficients,  the  STRIP  data  was 
reapplied  to  the  function  SECOND_SEARCH . M .  SECOND_SEARCH . M 
did  the  seune  thing  as  FIRST_SEARCH.M  but,  the  initial  guesses 
were  produced  by  the  coefficients  and  the  function  only  went 
through  one  minimization  routine  that  varied  all  of  the 
element  values.  Figure  17  depicts  the  flow  through  the  entire 
process . 
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Figur*  17.  The  entire  process  of  developing  the  coefficient 
from  the  raw  data  is  depicted  above. 
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EESOF  MODEL 

!  USER:  Michael  Linzey 

!  DATE:  3  MAR  91 

!  CIRCUIT:  model 

!  COMMENT:  This  is  a  model  that  represents  the  inductive  strip 

!  in  f inline.  The  element  values  must  be  entered  for  a 

!  specific  geometry. 

DIM 

FREQ  GHZ 
RES  OH 
IND  NH 
CAP  PF 
LNG  MIL 
TIME  PS 
COND  /OH 
ANG  DEG 

VAR 

!  GEOMETRY  SPECIFIC 
a=900 
b=400 
W=200 
d=0 
T=10 

!  ELEMENTS 

l_d»12.6320 
c_d=.0064 
turnO=l . 6020 
turnl=3 . 1142 
turn2=-2 . 5861 
turn3=0 .1429 

EQN 

PI=3. 14159 
A_a=a/2 

.•MODEL 

balance=l 

ba=b/a 

wb=W/b 

dal=d*2/a 

da=d/a 

•  EQUIVELENT  VALUES  FOR  THE  LOAD  HALF  WAVEGUIDE 
er_dl*. 999851+dal* . 0149217-dal**2* . 3391197 
er_d=er_dl+dal**3*10 . 61754-dal**4*3 . 664746+dal**5* . 135657 
A_dl=l+dal* . 001314-dal**2* . 0302559+dal**3* . 342420 
A  d=(A  dl-dal**4*2 . 161754+dal**5*l .7460658) *a/2 
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c=2. 998*10**8*1000/. 0254 
mu=4*pi*10** (-7) / (1000/. 0254) 
f=FREQ*10**9 
k0=2*pi*f /c 

!  beta  presumed  imaginary 
beta=sqrt (abs (er_d*k0**2- (pi/A_d) **2) ) 
kd=sqrt (2 . 22*k0**2+beta**2) 
ka=sqrt (k0**2+beta**2) 

CS=cos (d*kd) 

SN=sin (d*kd) 

zov_n=2*b*ka*kd**3* (2*pi*f ) *mu*CS**2 
zovl='ka**3*  (SN*CS-d*kd)  *cos  (A_d*ka-d*ka)  **2 
zov2=kd**3*CS**2*sin (A_d*ka-d*ka) *cos (A_d*ka-d*ka) 
zov_d= (beta* (zovl+zov2-ka*kd**3* (A_d-d) *CS**2) ) 
z  ov=z  o v_n / z  ov_d 

B_d=-zov*sqrt (abs (4*A_d**2*er_d*f**2-c**2) ) / (480*f*pi) 

!  CALCULATES  THE  FINLINE  EQUIVALENT  DIMENSIONS 
cl=-4 . 9723*ba**2+4 .7413*ba-0.7651 

Aeq= (2- (1- (ba+ . 45) * (1-wb) **2) **0 .5+cl* (1-wb) **26) *a 

c9=-20 . 16*da**2+6 . 42*da+ . 6494 
erl=c9* (1-wb** (1-exp (-10*da) ) ) 
er='erl+wb+2 . 604*da+  (1-da)  **6*  (1-wb) 

c2= (-115. 79*da**2+27. 87*da-. 4933) *ba+87 . 52*da**2-22 . 49*da- . 1932 
c3=0 . 29+0 . 0773*exp (l-40*da) 

c4*(20. 1154*da**2-3.729*da-0.0611) *ba+ (-26 . 1788*da**2+5 . 537*da+l . 0376) 
c5*-13.5217*da**2+2 .4017*da+0.0411 

Beq_avel=c2* (1-wb** (2*ba*c3) ) +c4+c5* (1-abs (ba-wb) **2) **4 

Bec^ave®  (Be<^_avel-0 . 025*  (1-abs  ( .  925-wb)  **2)  **16)  *b 

c6= (-76.251* (da) **2+17 . 23*da- . 1578) *ba+lll . 2*da**2-20 . 84*da- . 2936 

c7= (64 . 82*da**2-14 . 77*da- . 3029) *ba-107 . l*da**2+22 . 85*da- . 2936 

c8= (9 . 696*da**2-l . 449*da- . 1431) *ba-12 . 13*da**2+l . 39*da+ . 1195 

m=c6*wb**2+c7*wb+c8 

fc=(29980/2.54) / (2*a) 

Beq=m* (FREQ/ fc-1 .56) *b+Beq_ave 

! IMPEDANCE  TRANSFORMER 
Lambda= (29980/2.54) /FREQ 

LpovLl=l/ (1- (Lambda/ (2*Aeq*er**0 . 5) ) **2) **0.5 
LpovJ.2*l/  (1-  (Lambda/  (2*A) )  **2)  **0.5 
Z1=120*PI* (2*Beq/Aeq) *LpovLl/er**0 . 5 
f_n= (FREQ-fc) /fc 

Z_d=1000/«d3s (turn0+turnl*f_n+turn2*f_n**2+turn3*f_n**3) 
Z_a=Z_d/balance 

X_a=(Zl/Z_a)  **0.5 
X  d=(Zl/Z  d) **0.5 
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CKT 


IND  1  0  L''L_d 

CAP  1  0  C''C_d 

RWG  1  2  A''A_a  I.^'T  ER=1  RH0=1 

IND  2  0  L''L_d 

CAP  2  0  C''C_d 

DEF2P  1  2  Air 

IND  1  0  L''L_d 

CAP  1  0  C''C_d 

RWG  1  2  A''A_d  B^B_d  L'^T  ER''er_d  RH0=1 

IND  2  0  L''L_d 

CAP  2  0  C''C_d 

DEF2P  1  2  Dielec 

XFER  1200  N^X_a 

XFER  1  3  0  0  N''X_d 

Air  2  4 

Dielec  3  5 

XFER  6  5  0  0  N''X_d 

XFER  6  4  0  0  N''X_a 

DEF2P  1  6  STRIP 

RWGT  1  A^'Aeq  B'^Beq  ER'^er  RHO=l 

DEFIP  1  WEDGE 

!  SPECTRAL  DOMAIN  DATA 

S2PA  120  t50d30w5.S2P 

DEF2P  1  2  SPEC_DOM 

SPEC_DOM  1  2 
DEF2P  1  2  SPEC  M 


TERM 

PROC 


STRIP  WEDGE  WEDGE 


OUT 

STRIP  MAG(S11]  GRl 

STRIP  ANG[S11]  GRIA 

SPEC_DOM  MAG[S11]  GRl 
SPEC_DOM  ANG[S11]  GRIA 

FREQ 

SWEEP  8  12  1 


GRID 


RANGE 

8 

12 

1 

GRl 

0 

1 

.  1 

GRIA 

80 

180 
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The  following  listing  illustrates  the  form  of  the  spectral  domain  input 
data  file  for  the  iductive  strip  model.  This  is  the  listing  called  by 
the  model  on  the  previous  page. 


!  Spectral  Domain  data  file  for  the  following  geometry 
!  T=50  d=30  W=200 


!Freq 

ISllI 

<S11 

IS12I 

<S12 

IS21I 

<S21 

IS22I 

<S22 

8.0 

0.8992 

146.6191 

0.1914 

56.6191 

0.1914 

56.6191 

0 

.  8992 

146.6191 

9.0 

0.8528 

139.1836 

0.2727 

49.1836 

0.2727 

49.1836 

0 

.8528 

139.1836 

10.0 

0.8056 

132.0117 

0.3510 

42.0117 

0.3510 

42.0117 

0 

.8056 

132.0117 

11.0 

0.7568 

126.7910 

0.4272 

36.7910 

0.4272 

36.7910 

0 

.7568 

126.7910 

12.0 

0.7093 

121.3066 

0.4968 

31.3066 

0.4968 

31.3066 

0. 

7093 

121.3066 

13.0 

0.6610 

115.7168 

0.5630 

25.7168 

0.5630 

25.7168 

0 

.6610 

115.7168 
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The  following  is  a  listing  of  MATLAB  functions  used  in  element  calculations 

function  [f f mag, ang] =structure (fl, a,b, d, W, t, l_d, c_d, tO, tl, t2, t3) 

%  This  function  replicates  the  performance  of  TSTONE  software  for  the 
%  configuration  under  consideration. 

%  calculates  constants 
a_a=a/2;  er_a«l;  wb=W/b; 
fc=2 . 998e8*le-9*1000/ . 0254/ (2*a) 

%  iterates  through  frequencies 

for  n=l : length (fl) 

freq=fl (n) /  %  get  first  freq 

f=freq*le9; 

%  get  equivalent  loaded  waveguide 
[a_d, b_d, er_d] =half_guide (a, b, d, f ) ; 

yO_a=yov(a_a,b,er_a, f) ;  %  calculate  the  susceptance  of 

y0_d=yov (a_d,b_d, er_d, f ) ;  %  the  two  guides 

gamma_a=real ( j*beta (f , a_a, er_a) ) ;  %  calculate  the  prob.  const. 
gamma_d=real ( j*beta (f , a_d, er_d) ) ; 

fn= (freq-fc) /fc;  %  calculate  the  susptance  used  in 

y0_f=0 . 001*abs  (t0+tl*fn+t2*fn''2+t3*fn''3)  ;  %  turns  ratio 

y0_a=y0_a/y0_f ;  %  normalizes  the  impedance 

y  0__d=y  0_d/y  0_f ; 

%  pi  equivilant  for  waveguide 
ya=yO_a/sinh  (gaunma_a*t)  +y0_d/sinh  (gamma_d*t)  ; 
yb*yO_a*tanh (gamma_a*t/2) +y0_d*tanh (gamma_d*t/2) ; 

%  tank  cuircuits 
tank_a= ( j*2*pi*f *c_d+l/ ( j*2*pi*f*l_d) )/y0_f; 
tank_d=( j*2*pi*f*c_d+l/ ( j*2*pi*f *l_d) ) /y0_f; 

yb=yb+tank_a+tank_d;  %  form  the  complete  pi  oquivilant 

y_eq=l/ (1/ (1+yb) +l/ya) +yb;  %  find  the  eqpiivilent  suseptance 
sll= (l-y_eq) / (l+y_eq) ;  %  calculate  sll 

g  (n,  ; )  =  [f  *le-9  2Qds(s11)  angle  (sll)  *180/pi]  ; 
end 

f=g  ( ; , 1) ; 
mag=g ( : , 2 ) ; 
ang=g ( : , 3 ) ; 
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function  structure_l (f 1 , a, b, d, W, t) 

%  This  function  coupled  with  STRUCTURE_2  produce  the  saune  results  as 
%  STRUCTURE.  STRUCTURE_1  and  STRUCTURE_2  together  repoduce  the 
%  performance  of  TSTONE  for  the  device  under  concideration .  STRUCTURE_1 
%  computes  the  geometry  specific  parts  of  STRUCTURE  to  prevent 
%  for  each  frequency  iteration. 


the  freq[uencies 
(in  GHz) 


%  constants 

a_a=a/2;  n=0;  er_a=l;  wb=W/b; 
fc=le-9* (2 . 998e8*1000/ . 0254) / (2*a) ; 

%  iterate  through 

for  if req=l : length (fl) 

freq=fl (ifreq) ;  %  select  first  freq 

f=freq*le9; 

[a_d, b_d, er_d] =half_guide (a, b, d, f ) ; 
y0_a (if req) =yov (a_a, b, er_a, f ) ; 
y0_d  (ifreq)  =yov  (a_d,b_d,  er_d,  f )  ; 
gamma_a=real  ( j*beta  (f ,  a_a,  er_a)  )  ; 
gamma_d=real ( j*beta (f , a_d,  er_d) )  ; 
gamma_check (ifreq) =gamma_d; 

% 

sinh_a  (if  req)  =sinh  (gaunma_a*t )  ; 
sinh_d  (ifreq)  =sinh  (gaunma_d*t)  ; 
tanh_a (if req) =tanh (gamma_a*t/2) ; 
tanh_d (ifreq) =tanh (gamma_d*t/2) ; 
end 
f=fl; 

%  saves  calculations  to  CONSTANTS 
save  constants  f  yO_a  yO  d  gamma_check  fc  sinh  a  sinh  d 


%  get  the  equivilent  loaded  guide 
%  calculate  the  impedances  for 
the  two  guides 
calculate  the  probigation 
constants 

used  to  check  if  guide  starts 


% 

% 

% 

% 


to  probigate 

%  calculates 


need  trig  fucntions 


tanh  a  tanh  d 


function  [ f , mag, ang] =structure_2 (l_d, c_d, t0,tl,t2,t3) 

%  This  function  completes  the  calcuations  started  in  STRUCTURE_1 ;  includes 
%  the  effects  of  the  elements. 

load  constants  %  gets  constants  from  STRUCTURE_1 

l_d=abs (l_d) ;  c_d=abs (c_d) ;  n=l;  fl=f;  %  takes  the  abs  of  L  and  C  to 

%  prevent  the  use  of  negative  #'s 
for  freqi=l : length (f)  %  iterates  through  frequency 

if  gamma_check (n) ==0,  %  checks  to  see  if  the  guide  is 

%  probigating,  skips  if  yes 

else 

freq=f 1 (freqi) ;  f=freq*le9;  w=2*pi*f;  %  variables 

%  calculates  transformed  impedance 

f_n= (freq-fc) /fc; 

y0_f=0 . 001*abs  (t0+tl*f_n+t2*f_n''2+t3*f_n''3)  ; 

%  Normalizes  the  suseptance  of 
%  guides 

y0_a (n) =y0_a (n) /y0_f ; 
y0_d (n) =y0_d (n) /y0_f ; 
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%  pi-equivelant  of  guides 
ya=yO__a  (n)  /sinh_a  (n)  +yO_d  (n)  /sinh_d  (n) ; 
yb=yO_a (n) *tanh_a (n) +yO_d (n) *tanh_d (n) ; 

%  suseptance  of  the  tank  circuits 
tank_a= ( j*w*c_d+l/ ( j*w*l_d) ) /yO_f; 
tank_d= ( j*w*c_d+l/ ( j*w*l_d) ) /yO_f ; 

%  total  circuit 

yb*yb+ 1  ank_a + 1  ank_d ; 

%  eqivilent  suseptance 
y_eq=l/ (1/ (1+yb) +l/ya) +yb; 

%  calculate  sll 
sll=(l-y_eq) / (l+y_eq) ; 

g (n, : ) = [f *le-9  abs(sll)  angle (sll) *180/pi] ; 
n=n+l; 
end 
end 

f=g  ( : , 1) ; 
mag=g ( : , 2 ) ; 
ang=g( : , 3) ; 


The  following  functions  support  the  'structure'  functions: 
function  b=beta (f , a, er) 

%calcualtes  beta  for  a  given  frequency  for  the  TEIO  mode.  Assumes  that 
%a  is  measured  mils. 

c=2 . 998e8*1000/ . 0254/ 
k=2*pi*f*  (er)  ''.5/c/ 
kc=pi/a; 

b=(k^2-kc''2)  ''.5; 
end 


function  [a  d, b_d, er_d] =half_guide (a,b, d, f ) 

%This  functXon  computes  the  effective  dimentions  of  a  the  below  cutoff 
%waveguide  section  that  contains  dielectric.  The  model  assumes  that 
%the  dielectric  is  er=2.22. 

%  constants 

c=2.998e8*1000/.0254;  mu=pi*4e-7/ (1000/ . 0254) ;  d=d*2/a; 

%  calc  er 

er_d=.  999851+d*  .  0149217-d‘'2*  .  3391197  hd''3*10 . 61754 9-d''4*3 . 66474 e+d-'S*  .  135657 

r 

a_d=a/2*  (1+d*  .  001314-d''2*  .  0302559+d^3*  .  342420-d''4 *2 . 161754+d^5*l .  74  60658)  / 

%  calculate  Zov 

k0=2*pi*f /c; 

beta=sqrt  (er_d*k0^2- (pi/a_d)  ''2)  ; 
kd=sqrt  (2 . 22*k0''2-beta''2)  ; 
ka=sqrt  (k0''2-beta''2)  ; 
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CS=cos (d*kd)  ; 

SN=sin (d*kd)  ; 

zov_n=2*b*ka*kd'‘3*  (2*pi*f)  *mu*CS'^2; 
zovl=ka''3*  (SN*CS-d*kd)  *cos  (a_d*ka-d*ka)  ^2; 
zov2=kd''3*CS^2*sin  (a_d*ka-d*ka)  *cos  (a_d*ka-d*ka)  ; 
zov_d=  (beta*  (zovl+zov2-ka*kd''3*  (a_d-d)  *CS'^2) ) ; 
z  o v=  z  o v_n / z o v_d ; 

b  d=real  (-zov*sqrt  (4*a_d^2*er_d*f  ^2-c'^2)  /  (480*f  *pi)  )  ; 


function  y=yov (a, b, er, f ) 

%  calculates  voltage  power  suseptance 
fc=  (1000/  .0254)  *2.998e8/  (2*a*er'^0 . 5)  ; 
lambda* (2 . 998e8*1000/ . 0254) /f; 
z=(376.7/er'"0.5)  *  (2*b/a)  /  (l-(fc/f)  "'2)  '^0.5; 
y= (1/z)  ; 
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The  following  functions  are  connected  with  the  search  routines: 

function  error=calc_error (mag, s_mag, phase, s_phase, f , x) 

%  This  function  caclulates  the  error  between  the  two  input  magnitudes 
%  and  phases . 

s_mag=s_mag (1 : length (mag) ) ;  %  sets  the  length  of  the  vectors  to  be 

s_phase=s_phase  (1 :  length  (mag) )  ;  %  the  seune,  needed  in  the  case  of 

%  propagation  less  than  13GHz 

m=  (mag-s_mag)  ; 
p= (phase-s_phase) ; 

for  i=l : length (m)  %  finds  the  relative  error  at  each  point 

m(i)=m(i)  /s_mag(i)  ; 
p(i)=p(i)  /s_phase(i)  ; 
end 

global_m=max (abs (m) ) *100; 
global_p=max (abs (p) ) *100; 
error_mag=m' *m; 
error_jph=p'  *p; 

%  calculates  the  error  returned 
error=error_mag+error_ph-f  10*  (error_mag-error_ph)  ''2; 

%  viewing  options  when  running  in 
%  foreground 
%[global_m  global_p  x' ] 

%plot (f , s_mag, ' r' , f , s_phase/180, ' g' , f ,mag, ' *r' , f, phase/180, ' *g' ) 


function  error=model (x) 

%  This  function  takes  input  from  the  matlab  function  minimization 
%  routine  passes  values  to  STRUCTURE_2 .  STRUCTURE_2  calculates  the 
%  response  and  return  the  values  of  the  Sll  scattering  parameters. 

%  CALC_ERROR  then  compares  the  two  responses  and  calculates  an  error 
%  which  is  returned  to  the  function  minimization  routine,  for  use  in  the 
%  production  of  the  next  guess. 

load  sll_data_m 

l_d=x (1) *le-9;  %  element  values 

c_d=x (2) *le-12; 

turn0=x (3) ; 

turnl=x (4)  ; 

turn2=x (5)  ; 

turn3=x ( 6)  ; 

%  finds  response  for  given  element  values 
[f  mag  phase] =structure_2 (l_d, c_d, turnO, turnl , turn2 , turn3) ; 

%  calculates  error 

error=calc_error  (mag,  s_mag,  phase,  s__phase,  f ,  x)  ; 


function  error=model_2 (x) 

%  This  function  takes  input  from  the  matlab  function  minimization 
%  routine  passes  values  to  STRUCTURE  2.  STRUCTURE  2  calculates  the 
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%  response  and  return  the  values  of  the  Sll  scattering  pareimeters . 

%  CALC_ERROR  then  compares  the  two  responses  and  calculates  an  error 
%  which  is  returned  to  the  function  minimization  routine,  for  use  in  the 
%  production  of  the  next  guess . 
load  sll_data_m 
load  g_hold 

l_d=g_hold (1) *le-9;  %  element  values 

c_d=g_hold(2) *le-12; 

turnO=x (1) ; 

turnl=x (2) ; 

turn2=x (3)  ; 

turn3=x (4)  ; 

%  finds  response  for  given  element  values 
[f  mag  phase ] =structure_2 (l_d, c_d,turnO,turnl, turn2,turn3) ; 

%  calculates  error 

error=calc_error (mag, s_mag, phase, s_phase, f , x) ; 


function  g=f irst_search (g_init, tol, begin) 

%  This  function  finds  and  stores  element  values  that  match  the  scattering 
%  pareuneters  found  in  the  file  S11_DATA.  The  data  in  S11_DATA  is  in 
%  column  form  as  follows: 

%  [  freq,  sll_maq,  Sll_phase,  t,  d,  W] 

%  The  data  is  first  broken  in  to  blocks  that  represent  each  geometry.  The 
%  data  is  then  saved  in  Sll_data_m  for  recall  by  the  function  minimization 
%  routines.  STRUCTURE_1  then  calculates  the  geometry  dependant  parameters. 
%  Initail  starting  points  for  element  values  are  eather  given  by  the 
%  argument  or  from  the  previous  iteration.  Function  minimiztion  routines 
%  are  then  called  to  find  values  of  the  elements  that  match  the  S11_DATA. 

%  The  variables  global_m  and  global_p  must  be  declared  to  be  global . 

%  initial  values 

g=g_init;  g_select=l;  d_hold=0;  a=900;  b=400;  n=0; 
load  sll_data 

%  block  data  by  geometry,  using  the 
%  starting  frequency 
for  i=l : length (sll_data) 
if  7  ==  sll_data (i, 1)  , 
n=n+l ; 

block (n)=i;  %  record  begining  of  each  block 

end 
end 

block  (n+1)  =length  (sll_data) +1;  %  mark  the  end  of  the  last  block 

for  i=begin : length (block) -1  %  iterate  through  all  the  data 

i_end=block (i+1) -1;  %  calculate  the  end  of  the  block 

f=sll_data (block (i) +1 : i_end, 1 ) ;  %  block  out  the  data 

s_mag=sll_data (block (i) +1 : i_end, 2) ; 
s_phase=sll_data (block (i) +1 : i_end, 3) ; 
t=sll_data (block (i)  ,4)  ; 
d=sll  data (block (i)  ,  5)  ; 
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%  test  to  see  if  magnitude  of 
%  sll  is  increasing 


W=sll_data (block (i) , 6) ; 
i_end*length (s_mag) ; 
for  k*length(s_mag)-l:-l;2 
if  s_mag (k) <s_mag (k+1) 
i_end=k; 
end 
end 

f =f ( 1 : i_end) ;  %  eliminate  the  bad  points  noted 

s_mag*s_mag  ( 1 :  i_end)  ;  %  eUbove 

s_j>hase=s_phase  (1 ;  i_end)  ; 

(d  t  W]  %  display  block  being  worked  on 

[f  s_mag  s_phase] 
save  sll_data_m 

structure_l (f , a,b, d, W, t) ;  %  calculate  geometry  dependant 

%  values 

g_hold*g  %  save  guess  of  element  values 

save  g_hold  g_hold 

%  pass  turns  ratio  into  function 
%  minimization  routine 

g_hold=fmins  ('model_2' ,  [g  (3)  g(4)  g(5)  g(6)],0.10); 

%  re-form  element  vector 

g=tg(l)  g(2)  g_hold(l)  g_hold(2)  g_hold(3)  g_hold(4)] 
g_hold*g 

[global_m  global_p] 
save  g_hold  g_hold 

[g,count]=fmins ('model' ,g,tol) ;  %  pass  all  elements  into  function 

%  miniminzation  routine 

store (a, b, W, d, t , tol, g' , count/100)  %  save  the  element  values 

if  t<15, 

g_save=g  %  save  the  value  of  the  begining 

end  %  of  one  d  set  to  start  the  next 

if  t>950 
g=g_save 
end 
end 


function  g=second  search (tol, begin) 

%  Operation  of  thTs  function  is  the  same  as  the  FIRST_SEARCH  except 
%  in  the  way  the  initial  guess  of  the  elements  is  found.  As  note 
%  below. 

n*0;  a=900;  b*=400; 
load  sll__data 
for  i=l : length (sll_data) 
if  7  ==  sll_data (i, 1)  , 
n=n+l ; 
block (n) =i; 
end 
end 
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block  <n+l) “length (sll_data) +1; 
for  i“begin; length (block) -1 
i_end“block (i+1) -1; 
f=3ll_data (block (i) +1 ; i_end, 1) ; 
s_inag=sll_data  (block  (i)  +1 :  i_end,  2)  ; 
s_phase=sll_data (block (i) +1 ; i_end, 3) ; 
t=sll_data (block {i) ,4) ; 
d“Sll_data (block (i) ,5) ; 

W=sll_data (block (i) , 6) ; 
i_end=length (s_mag) ; 
for  k=length(s_mag) -1;-1:2 
i  f  s_inag  ( k )  <s_mag  ( k+ 1 ) 
i_end“k ; 
end 
end 

f=f (1 : i_end) ; 
s_mag=s_mag ( 1 : i_end) ; 
s_phase=s_phase (1 ; i_end) ; 

[d  t  W] 

(f  s_mag  s_phase] 
save  sll_data_m 
structure_l (f , a, b, d, W, t ) ; 

%  gets  the  elements  as  a  function 
%  of  parameters 

[L,C, to, tl,t2, t3]=get_elements (a,b,d,W, t) ; 
g=[L  C  to  tl  t2  t3] 

[g, count ] “fmins ('model' ,g,tol) ; 
store2  (a,b,  d,  t ,  tol,  g' ,  count/100) 
if  t<15, 
g_save*g 
end 

if  t>950 
g=g_save 
end 
end 
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The  following  functions  relate  to  coefficient  reduction  and  element 
generation . 


function  coef_d=find_coef (s) 

%  This  function  find  the  coefficients  of  the  polynomials  using  MATLAB's 
%  mean  square  function. 

Nt=2;  Nd=2;  Nb='2; 
s (10 : Ilf : ) =abs (s (10 : 11 , : ) ) ; 
s=[s  zeros (16, 1) ] / 
for  n=l ; length (s (5, : ) ) -1 
if  s (5, n) >s (5, n+1) 

Ne=n; 

for  i=10:15  %  finds  the  coef  of  the  elements 

%  with  respect  to  T/a 

coef_t= [coef_t ;  polyfit (s (5,Nb;Ne) /s (1, 1) , s (i,Nb:Ne) ,Nt) '  ]  ; 
end 

coef_l=(coef  1  coef_t]  %  puts  the  results  in  a  column 

d=[d  s (4,Nb)7s (1, 1) ]  %  makes  d/a  list  to  reduce  the 

%  above  list 

coef_t= [ ] ; 

Nb=Ne+l; 

end 

end 

[row, col] =size (coef_l) ; 

for  n=l:row  %  finds  the  coef  of  the  coef 

%  generated  above  wrt  d/a 
coef_d«  [coef__d;  polyf  it  (d,  coef_l  (n,  : ) ,  Nd)  '  J 
end 


function  coef_wb=f ind_wb (wb, coef_in) 

%  Finds  the  coefficients  with  respect  to  W/b 
[r, c] =size (coef_in) / 

Nw=2  ; 
for  i=l:r 

coef_wb= [coef_wb;  polyf it (wb, coef_in (i, ; ) ,Nw) ] ; 
end 
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function  [L, C, tO, tl, t2, t3] =get_elements (a,b, d, W, t) 

%  This  function  used  the  values  of  the  coefficients  in  the  file  COEF 
%  to  calculate  the  element  values . 

load  coef 

Nd=2;  %  the  order  of  the  d/a  polynomial 

%  needed  for  initial  decomposition 
coef_d=get_wb (Nt , coef , W/b) ;  %  accounts  for  the  W/b  variation 

(r, c] =size (coef_d) ; 

N=r/6;  da=d/a; 

for  i=l:N  %  iterates  throught  the  coef;  accounts  for 

%  the  d/a  variation 

coef_L= [  coef_L  polyval (coef_d (i, : ) , da)  ] ; 
coef_C= [  coef_C  polyval (coef_d (N+i, : ) , da)  ] ; 
coef_tO=[  coef_tO  polyval (coef_d (2*N+i, : ) ,da)  ]  ; 
coef_tl= [  coef_tl  polyval (coef_d (3*N+i, : ) , da)  ] ; 
coef_t2=[  coef_t2  polyval (coef_d (4 *N+i, :), da)  ]; 
coef_t3=[  coef_t3  polyval (coef_d (5*N+i, ;), da)  ]; 
end 

%  computes  the  value  of  the  elements 
L=polyval (coef_L,t/a) ; 

C=polyval (coef_C,t/a)  ; 
tO=polyval (coef_tO, t/a) ; 
tl=polyval (coef_tl , t/a) ; 
t2=polyval (coef_t2, t/a) ; 
t3=polyval (coef  t3,t/a) / 


function  coef_d=get  wb (Nb, coef_wb, wb) 

%  Account  for  the  w7b  variation 
for  i=l ; length (coef_wb) 

cO (i) =polyval (coef_wb (i, ; ) , wb) ;  %  calculates  the  W/b  variation 

end 

for  i*l ;Nd+l : length (cO)  %  resizes  the  matrix  for  next 

cl=[  cl;  cO(i:i+Nd)  ];  %  operation 

end 

coef  d=cl; 


59 


APPENDIX  D.  COEFFICIENT  REDUCTION  OUTPUT 


Following  this  page  are  three  sample  outputs  that 
illustrate  the  difficulties  encountered  during  the  development 
of  the  coefficients  which  were  to  represent  the  element 
values.  The  first  listing  shows  the  results  of  a  function 
minimization  where  the  inductance  was  allowed  to  vary  and  the 
capacitance  and  the  four  terms  of  the  turns  ratio  were 
calculated  by  the  coefficients.  The  second  listing  shows  the 
percent  error  in  the  desired  value  of  the  inductance  and  the 
value  produced  from  the  least  mean  square  generating  programs. 
The  third  listing  shows  the  percentage  error  in  phase  and 
magnitude  with  the  STRIP  data  when  the  coefficients  calculated 
above  were  used  to  generate  element  values  in  the  model . 
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The  following  listing  resulted  from  generating  element  values  for 
all  of  the  elements  except  the  inductor  using  coefficients  and  then 
using  a  function  minimization  routine  to  generate  the  inductor 
value.  The  columns  are  arranged  as  follows; 
a 
b 
W 
d 
T 

tolerance 

(search  count) *0.01 
%  phase  error 
%  mag  error 
L  inductor  in  nH 

C  capacitor  in  pF 

to  0  order  coefficient  for  the  turns  ratio 
tl  1st  order  coefficient  for  the  turns  ratio 

t2  2nd  order  coefficient  for  the  turns  ratio 

t3  3rd  order  coefficient  for  the  turns  ratio 


Columns  1  through  7 


900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

400 . 0000 

400 . 0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

0 

0 

0 

0 

0 

0 

0 

10.0000 

20.0000 

30.0000 

40.0000 

50.0000 

60.0000 

70.0000 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0900 

0.0900 

0.0900 

0.0900 

0.0900 

0.0900 

0.0900 

2.6978 

1.8785 

0.9913 

0.9183 

0.8394 

1.4123 

2.0688 

0.4667 

0.5811 

0.8052 

1.0324 

1.2725 

1.6186 

1.9703 

17 . 6203 

17 . 0645 

16.7005 

16.4661 

16.3048 

16.1440 

16.0190 

0.0047 

0.0049 

0.0052 

0.0054 

0.0057 

0.0060 

0.0063 

1.5539 

1.4595 

1 . 3682 

1.2798 

1.1946 

1.1123 

1.0331 

2.5056 

2.9212 

3.3246 

3.7159 

4.0949 

4.4617 

4 . 8164 

1.5626 

0.5833 

-0.3683 

-1.2922 

-2.1885 

-3.0571 

-3.8981 

-3.2824 

-2.6031 

-2.1004 

-1.5344 

-0.9851 

-0.4524 

0.0636 

0 

0 

0 

0 

0 

0 

0 
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Coliimns  8  through  14 


900.00C0 
400 . 0000 
400.0000 
0 

80.0000 
0.0100 
0.0900 
2.5424 
2 . 3146 
15.8997 
0.0065 
0 . 9570 
5 . 1589 
-4.7114 
0 . 5630 
0 

Columns 

900 . 0000 
400 . 0000 
400.0000 
10 . 0000 
10 . 0000 
0.0100 
0.0900 
2.7920 
1.0093 
15.4365 
0 . 0047 
1.4917 
3.2785 
-0.6303 
-1 . 8272 
0 


900.0000 

400.0000 

400.0000 

0 

90.0000 
0.0100 
0.0900 
1 . 1586 
0.3406 
16.1901 
0 . 0068 
0.8838 
5.4892 
-5.4970 
1.0457 
0 

through 

900 . 0000 
400.0000 
400.0000 
10.0000 
20.0000 
0.0100 
0.0900 
1 . 9233 
1 . 1187 
14 . 7504 
0 . 0048 
1.4271 
3.4819 
-1.1668 
-1 . 4472 
0 


900.0000 

400.0000 

400.0000 

0 

100.0000 

0.0100 

0.0900 

1.1907 

0.6124 

16.1876 

0.0071 

0.8137 

5.8073 

-6.2549 

1.5117 

0 


900.0000 
400.0000 
400.0000 
10.0000 
30.0000 
0.0100 
0.0900 
1.2957 
1.2459 
14.2990 
0.0049 
1.3647 
3.6807 
-1 . 6912 
-1 . 0769 
0 


900.0000 

400.0000 

400.0000 

0 

200  0000 
0.0100 
0.1000 
0.6524 
0.9522 
16.5842 
0.0101 
0.2799 
8.3180 
-12.3128 
5.2560 
0 


900.0000 

400.0000 

400.0000 

10.0000 

40.0000 

0.0100 

0.0900 

0.5927 

1.3433 

13.9707 

0.0050 

1.3045 

3.8750 

-2.2036 

-0.7162 

0 


900.0000 

400.0000 

400.0000 

0 

300.0000 

0.0100 

0.1000 

0.9807 

0.8626 

15.8695 

0.0136 

0.0499 

9.6100 

-15.6042 

7.3344 

0 


900.0000 

400.0000 

400.0000 

10.0000 

50.0000 

0.0100 

0.0900 

0.5989 

1.4925 

13.7161 

0.0051 

1.2465 

4.0648 

-2.7039 

-0.3651 

0 


900.0000 

400.0000 

400.0000 

0 

400.0000 
0.0100 
0.0800 
3.7443 
0.0775 
13.2978 
0.0174 
0 . 1236 
9.6833 
-16.1289 
7.7469 
0 


900.0000 

400.0000 

400.0000 

10.0000 

60.0000 

0.0100 

0.0900 

0.5894 

1.6206 

13.5044 

0.0052 

1.1907 

4.2501 

-3.1922 

-0.0237 

0 


900.0000 

400.0000 

400.0000 

0 

500.0000 
0.0100 
0 . 0800 
2.2694 
0.2237 
9.6489 
0.0216 
0.5011 
8.5378 
-13.8871 
6.4936 
0 


900 . 0000 
400.0000 
400.0000 
10.0000 
70.0000 
0.0100 
0.0900 
0.7307 
1.7332 
13.3290 
0.0054 
1.1370 
4.4309 
-3.6683 
0.3081 
0 
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Columns  22  through  28 


900 . 0000 

900.0000 

900 . 0000 

900.0000 

900.0000 

900.0000 

900 . 0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

80.0000 

90.0000 

100.0000 

200.0000 

300.0000 

400.0000 

500.0000 

0 . 0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0900 

0.0900 

0.0900 

0.0900 

0.0800 

0.0800 

0.0900 

0.7781 

0.8612 

0.9399 

0.9337 

0.1064 

0.0487 

0.0699 

1.7554 

1.8820 

1.9826 

3.8107 

1.9144 

1.6289 

0.7899 

13.1603 

13.0097 

12 . 8710 

11.5923 

10.9954 

9.1540 

7.1753 

0.0055 

0.0057 

0.0058 

0.0079 

0.0109 

0.0148 

0.0196 

1 . 0855 

1.0362 

0.9890 

0.6374 

0.5035 

0.5874 

0.8890 

4.6071 

4.7788 

4.9460 

6.3695 

7.3412 

7.8610 

7.9289 

-4 . 1324 

-4.5845 

-5.0244 

-8.7601 

-11.2889 

-12 . 6106 

-12.7254 

0 . 6303 

0 . 9428 

1.2456 

3.7443 

5.2792 

5.8502 

5.4575 

0 

0 

0 

0 

0 

0 

0 

Columns  29 

1  through 

35 

900.0000 

900.0000 

900 . 0000 

900.0000 

900.0000 

900.0000 

900.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

10.0000 

20.0000 

30.0000 

40.0000 

50.0000 

60.0000 

70.0000 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0100 

0.0900 

0.0900 

0.0900 

0.0900 

0.0900 

0.0900 

0.0900 

1.0938 

1.0476 

1.1273 

0.9510 

0.9825 

0.9362 

1.0262 

0.4272 

0 . 6167 

0.7064 

0.7250 

0 . 6569 

0.5621 

0.3492 

15.4914 

14 . 9184 

14.5766 

14.3525 

14 .1980 

14.0834 

14.0077 

0.0049 

0 . 0048 

0.0047 

0.0046 

0 . 0045 

0.0045 

0 . 0045 

1.4964 

1 . 4717 

1.4480 

1.4251 

1.4032 

1.3821 

1.3620 

4.3038 

4.2170 

4.1360 

4.0607 

3.9911 

3.9273 

3.8692 

-3.7174 

-3.6376 

-3.5668 

-3.5051 

-3.4525 

-3.4091 

-3.3747 

0 . 1571 

0.2234 

0.2898 

0.3563 

0.4231 

0 . 4899 

0.5569 

0 

0 

0 

0 

0 

0 

0 
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Colunuis  36  through  42 


900.0000 
400.0000 
400.0000 
30.0000 
80.0000 
0.0100 
0.0900 
1.0220 
0.1549 
13.9314 
0.0045 
1.3428 
3.8168 
-3.3495 
0 . 6241 
0 

Columns 

900.0000 
400.0000 
400.0000 
50.0000 
10.0000 
0.0100 
0.0900 
0.4114 
0.6316 
15.5422 
0.0050 
1 . 6733 
4.6350 
-5.0731 
0.9066 
0 


900.0000 
400.0000 
400.0000 
30.0000 
90.0000 
0.0100 
0.0900 
0 . 9504 
-0.1042 
13.8818 
0.0045 
1 . 3245 
3.7701 
-3.3334 
0.6914 
0 

through 

900.0000 
400.0000 
400.0000 
50.0000 
20.0000 
0.0100 
0.0900 
0.5845 
0.9808 
14 . 9465 
0.0049 
1.6621 
4.4373 
-4.7357 
0.8258 
0 


900.0000 
400,0000 
400.0000 
30.0000 
100.0000 
0.0100 
0 . 0900 
1.0565 
-0.3019 
13.8258 
0.0045 
1.3071 
3.7292 
-3.3263 
0.7588 
0 


900.0000 
400.0000 
400.0000 
50.0000 
30.0000 
0.0100 
0.0900 
0.7866 
1 . 1108 
14,5972 
0.0049 
1 . 6512 
4.2488 
-4.4152 
0.7492 
0 


900.0000 

400.0000 

400.0000 

30.0000 

200.0000 

0.0100 

0.0800 

0.6479 

0.5094 

13.2444 

0.0055 

1.1833 

3.6353 

-3.7568 

1.4416 

0 


900.0000 

400.0000 

400.0000 

50.0000 

40.0000 

0.0100 

0.0900 

0.8306 

1.2071 

14.3616 

0.0048 

1.6406 

4.0694 

-4.1116 

0.6769 

0 


900.0000 

400.0000 

400.0000 

30.0000 

300.0000 

0.0100 

0.0800 

0.3448 

0.7133 

11.7661 

0.0081 

1.1507 

4.1144 

-5.0978 

2.1394 

0 


900.0000 

400.0000 

400.0000 

50.0000 

50.0000 

0.0100 

0.0900 

0.9956 

1.1559 

14.2099 

0.0047 

1.6302 

3.8993 

-3.8249 

0.6088 

0 


900.0000 

400.0000 

400.0000 

30.0000 

400.0000 

0.0100 

0.0800 

0.2139 

0.5648 

9.7564 

0.0123 

1.2092 

5.1667 

-7.3493 

2.8521 

0 


900.0000 

400.0000 

400.0000 

50.0000 

60.0000 

0.0100 

0.0900 

1.0635 

1.0530 

14.1129 

0.0047 

1.6200 

3.7384 

-3.5552 

0.5450 

0 


900.0000 

400.0000 

400.0000 

30.0000 

500.0000 

0.0100 

0.0800 

0.1880 

3.6981 

7.8159 

0.0181 

1.3589 

6.7920 

-10.5114 

3.5799 

0 


900.0000 
400.0000 
400.0000 
50.0000 
70.0000 
0.0100 
0.0900 
1.0886 
0.8981 
14.0423 
0.0047 
1 . 6101 
3.5867 
-3.3023 
0.4855 
0 
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Columns  50  through  56 


900.0000 
400 . 0000 
400.0000 
50.0000 
80.0000 
0.0100 
0.0900 
1 . 1351 
0.7281 
13.9922 
0.0047 
1.6004 
3.4442 
-3.0665 
0.4303 
0 

Columns  i 

900.0000 
400 . 0000 
400.0000 
70 . 0000 
10.0000 
0.0100 
0.0900 
0.9719 
0.6034 
15.3784 
0 . 0052 
2.0224 
4.2721 
-4.6973 
0.4213 
0 


900.0000 
400.0000 
400.0000 
50.0000 
90.0000 
0.0100 
0.0900 
1 . 1735 
0.4689 
13.9515 
0.0047 
1.5910 
3.3108 
-2 . 8475 
0.3793 
0 

through 

900.0000 
400.0000 
400.0000 
70.0000 
20.0000 
0.0100 
0.0900 
0.4303 
0.7000 
14.6776 
0.0053 
1.9983 
4 . 1426 
-4.4611 
0.3600 
0 


900.0000 

400.0000 

400.0000 

50.0000 

100.0000 

0.0100 

0.0900 

1.1892 

0.2374 

13.9251 

0.0048 

1.5819 

3.1867 

-2.6455 

0.3325 

0 
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900.0000 
400.0000 
400.0000 
70.0000 
30.0000 
0.0100 
0.0900 
0.5465 
0.7568 
1  '  .2433 
0.0054 
1 . 9745 
4.0190 
-4.2364 
0.3013 
0 


900.0000 

400.0000 

400.0000 

50.0000 

200.0000 

0.0100 

0.0800 

1.2297 

0.5446 

13.4752 

0.0060 

1.5040 

2.4512 

-1.5561 

0.0999 

0 


900.0000 
400.0000 
400.0000 
70.0000 
40.0000 
0.0100 
0.0900 
0.6036 
0.7414 
13.9576 
0.0056 
1 . 9509 
3.9013 
-4.0230 
0.2454 
0 


900.0000 
400.0000 
400.0000 
50.0000 
300.0000 
0.0100 
0.0800 
1.5884 
0.7871 
12.1279 
0.0090 
1.4512 
2 . 6354 
-2.1593 
0.2939 
0 


900.0000 

400.0000 

400.0000 

70.0000 

50.0000 

0.0100 

0.0900 

0.7790 

0.7860 

13.7379 

0.0057 

1.9275 

3.7895 

-3.8210 

0.1922 

0 


900.0000 
400.0000 
400.0000 
50.0000 
400.0000 
0.0100 
0.0800 
0 . 9332 
0.1534 
10.1444 
0.0136 
1.4235 
3.7394 
-4.4551 
0 . 9145 
0 


900.0000 

400.0000 

400.0000 

70.0000 

60.0000 

0.0100 

0.0900 

1.0208 

0.7618 

13.5720 

0.0059 

1.9043 

3.6835 

-3.6304 

0.1417 

0 


900.0000 
400.0000 
400.0000 
50.0000 
500.0000 
0.0100 
0.0800 
0 . 9515 
2.2812 
7.9998 
0.0198 
1.4208 
5.7630 
-8 . 4436 
1.9617 
0 


900.0000 
400.0000 
400.0000 
70.0000 
70.0000 
0.0100 
0.0900 
1.0358 
0.6843 
13.4432 
0.0061 
1 . 8812 
3.5834 
-3.4512 
0.0939 
0 
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Columns  64  through  70 


900.0000  900.0000  900.0000  900.0000  900.0000  900.0000  900.0000 

400.0000  400.0000  400.0000  400.0000  400.0000  400.0000  400.0000 

400.0000  400.0000  400.0000  400.0000  400.0000  400.0000  400.0000 

70.0000  70.0000  70.0000  70.0000  70.0000  70.0000  70.0000 

80.0000  90.0000  100.0000  200.0000  300.0000  400.0000  500.0000 

0.0100  0.0100  0.0100  0.0100  0.0100  0.0100  0.0100 

0.0900  0.0900  0.0900  0.0800  0.0800  0.0800  0.0800 

1.1151  1.2713  1.2846  1.2579  0.7262  0.5771  0.4570 

0.6889  0.6666  0.6126  0.6307  1.7171  2.5476  2.6323 

13.3299  13.2318  13.1585  12.5084  11.8911  10.5411  8.9369 

0.0063  0.0065  0.0067  0.0094  0.0134  0.0185  0.0248 

1.8584  1.8357  1.8133  1.5995  1.4051  1.2302  1.0749 

3.4892  3.4009  3.3185  2.8174  2.9043  3.5791  4.8419 

-3.2833  -3.1269  -2.9819  -2.1580  -2.4734  -3.9280  -6.5218 

0.0488  0.0064  -0.0332  -0.2808  -0.2573  0.0373  0.6029 

0  0  0  0  0  0  0 

Columns  71  through  77 

900.0000  900.0000  900.0000  900.0000  900.0000  900.0000  900.0000 

400.0000  400.0000  400.0000  400.0000  400.0000  400.0000  400.0000 

400.0000  400.0000  400.0000  400.0000  400.0000  400.0000  400.0000 

90.0000  90.0000  90.0000  90.0000  90.0000  90.0000  90.0000 

10.0000  20.0000  30.0000  40.0000  50.0000  60.0000  70.0000 

0.0100  0.0100  0.0100  0.0100  0.0100  0.0100  0.0100 

0.0900  0.0900  0.0900  0.0900  0.0900  0.0900  0.0800 

2.2465  0.8248  0.6130  1.0333  1.4566  1.7532  1.9895 

1.1545  1.2252  1.3701  1.4892  1.6555  1.9017  2.0748 

15.1502  14.2180  13.5937  13.1212  12.7360  12.4109  12.1343 

0.0054  0.0059  0.0064  0.0070  0.0075  0.0080  0.0086 

2.5436  2.4804  2.4179  2.3561  2.2952  2.2349  2.1754 

3.2151  3.3331  3.4468  3.5563  3.6615  3.7625  3.8594 

-2.5900  -2.8140  -3.0304  -3.2393  -3.4408  -3.6347  -3.8212 

-1.2988  -1.1741  -1.0538  -0.9381  -0.8268  -0.7201  -0.6179 

0  0  0  0  0  0  0 


66 


Colxonins  78  through  84 


900.0000 
400 . 0000 
400.0000 
90.0000 
80.0000 
0.0100 
0.0800 
3.5099 
1.4109 
12.5358 
0.0091 
2.1166 
3.9520 
-4 . 0001 
-0.5203 
0 


900.0000 

400.0000 

400.0000 

90.0000 

90.0000 

0.0100 

0.0800 

3.7978 

1.5650 

12.2834 

0.0097 

2.0586 

4.0403 

-4.1715 

-0.4271 

0 


900.0000 
400.0000 
400.0000 
90.0000 
100.0000 
0.0100 
0.0800 
3.9986 
1.7247 
12.0626 
0.0102 
2.0014 
4 . 1245 
-4.3355 
-0.3385 
0 


900.0000 

400.0000 

400.0000 

90.0000 

200.0000 

0.0100 

0.0800 

3.0324 

-0.6245 

10.5341 

0.0157 

1.4696 

4.7338 

-5.5627 

0.2995 

0 


900.0000 

400.0000 

400.0000 

90.0000 

300.0000 

0.0100 

0.0800 

2.1798 

0.1249 

9.7944 

0.0213 

1.0123 

4.9209 

-6.0402 

0.4859 

0 


900.0000 

400.0000 

400.0000 

90.0000 

400.0000 

0.0100 

0.0800 

0.6948 

2.7824 

10.0006 

0.0271 

0.6294 

4.6859 

-5.7679 

0.2205 

0 


900.0000 

400.0000 

400.0000 

90.0000 

500.0000 

0.0100 

0.0900 

0.3984 

6.8887 

9.9825 

0.0330 

0.3210 

4.0288 

-4.7460 

-0.4965 

0 
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The  following  is  a  partial  listing  of  the  output  of  a 
MATLAB  function  that  gives  the  percent  error  in  the  desired 
element  value  and  the  value  calculated  from  the  coefficients 
at  specific  geometries.  The  data  is  arranged  in  columns,  the 
first  five  places  represent,  a,b,W,d,T  and  the  last  five 
places  represent  the  element  values .  As  discussed  above  only 
the  inductor  shows  any  significant  deviation  from  the  desired 
values . 


element  error  = 


Columns  1  through  7 


900.0000 
400 . 0000 
400.0000 
0 

10.0000 
-17.3010 
-0.0000 
0 . 0000 
-0.0000 
0.0000 
0.0000 

Columns  8 

900.0000 

400.0000 

400.0000 

0 

80.0000 

-5.7242 

-0.0000 

-0.0000 

-0.0000 

0.0000 

-0.0000 


900.0000 

400.0000 

400.0000 

0 

20.0000 

-13.6287 

-0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 

through 

900.0000 

400.0000 

400.0000 

0 

90.0000 

-7.8888 

-0.0000 

-0.0000 

-0.0000 

0.0000 

-0.0000 


900.0000 

400.0000 

400.0000 

0 

30.0000 

-10.9892 

-0.0000 

0.0000 

-0.0000 

-0.0000 

0.0000 


900.0000 

400.0000 

400.0000 

0 

100.0000 

-8.5581 

-0.0000 

-0.0000 

-0.0000 

0.0000 

-0.0000 


900.0000 

400.0000 

400.0000 

0 

40.0000 

-9.1597 

-0.0000 

0 

-0.0000 

-0.0000 

0.0000 


900.0000 

400.0000 

400.0000 

0 

200.0000 

-28.6469 

0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 


900.0000 

400.0000 

400.0000 

0 

50.0000 

-7.9007 

-0.0000 

0 

-0.0000 

-0.0000 

0.0000 


900.0000 

400.0000 

400.0000 

0 

300.0000 

-65.4846 

0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 


900.0000 

400.0000 

400.0000 

0 

60.0000 

-6.8295 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

0.0000 


900.0000 

400.0000 

400.0000 

0 

400.0000 
-132.0765 
0.0000 
0.0000 
-0.0000 
0.0000 
0.0000 


900.0000 

400.0000 

400.0000 

0 

70.0000 

-6.1583 

-0.0000 

-0.0000 

-0.0000 

0 

-0.0000 


900.0000 

400.0000 

400.0000 

0 

500.0000 

-280.2833 

0.0000 

0.0000 

-0.0000 

0.0000 

0.0000 
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Columns  15  through  21 


900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900. 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400. 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400. 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10. 

10.0000 

20.0000 

30.0000 

40.0000 

50.0000 

60.0000 

70  . 

-5.0626 

-0.0575 

3.5428 

6.2640 

8.3605 

10.0122 

11. 

0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0. 

0.0000 

0.0000 

0 

-0.0000 

-0.0000 

-0.0000 

-0. 

-0.0000 

-0 . 0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0. 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

0.0000 

0.0000 

0. 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0. 

Columns  22 

1  through 

28 

900 . 0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900. 

400 . 0000 

400 . 0000 

400.0000 

400.0000 

400.0000 

400.0000 

400. 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400. 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10.0000 

10. 

80 . 0000 

90.0000 

100.0000 

200.0000 

300.0000 

400.0000 

500. 

12.2612 

12 . 9748 

13.4181 

6.2806 

-29.8648 

-91.5673 

-218 

-0.0000 

-0 . 0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0, 

-0.0000 

-0.0000 

-0.0000 

0.0000 

0.0000 

0.0000 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0 

Columns  29  through 

35 

900.0000 

900 . 0000 

900.0000 

900.0000 

900.0000 

900.0000 

900 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30 

10.0000 

20.0000 

30.0000 

40.0000 

50.0000 

60.0000 

70 

-4 . 9390 

-1.5317 

0.4749 

1.6842 

2.3735 

2.7318 

2 

0.0000 

0.0000 

0.0000 

0 

-0.0000 

-0.0000 

-0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0 

0000 

0000 

0000 

0000 

0000 

2380 

0000 

0000 

0000 

0000 

0000 


0000 

0000 

0000 

0000 

0000 

.4440 

.0000 

.0000 

.0000 

.0000 

.0000 


.0000 
.0000 
.0000 
.0000 
.0000 
.7550 
.0000 
.0000 
.  0000 
.0000 
.0000 
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Columns  36  through  42 


900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400 . 0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

30.0000 

80.0000 

90.0000 

100.0000 

200.0000 

300.0000 

400.0000 

500.0000 

2.7284 

2.4490 

2.1606 

-3.8652 

-10.0880 

-21.4418 

-48.9710 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0  0000 

-0.0000 

-0.0000 

-0.0000 

0 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 

-0.0000 
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The  following  is  a  listing  of  the  output  of  a  function 
that  iterates  through  the  STRIP  data  and  gives  the  resulting 
maximum  percent  error  for  that  geometry  with  the  model  using 
element  values  generated  using  coefficients. 

The  columns  are  arranged  in  the  following  way: 

a 

b 

w 

d 

T 

L  inductor  in  nH 

C  capacitor  in  pF 

to  0  order  coefficient  for  the  turns  ratio 
tl  1st  order  coefficient  for  the  turns  ratio 

t2  2nd  order  coefficient  for  the  turns  ratio 

t3  3rd  order  coefficient  for  the  turns  ratio 

%  error  phase 

%  error  magnitude 


Columns  1  through  7 


900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

0 

0 

0 

0 

0 

10.0000 

10.0000 

10.0000 

40.0000 

70.0000 

100.0000 

400.0000 

20.0000 

50.0000 

14.5784 

14.9579 

15.0325 

14.8023 

-4.2655 

14.7420 

14.8628 

0.0047 

0.0054 

0.0063 

0.0071 

0.0174 

0.0048 

0.0051 

1.5539 

1.2798 

1.0331 

0.8137 

0.1236 

1.4271 

1.2465 

2.5056 

3.7159 

4.8164 

5.8073 

9.6833 

3.4819 

4.0648 

1.5626 

-1.2922 

-3.8981 

-6.2549 

-16.1289 

-1.1668 

-2.7039 

-3.2824 

-1.5344 

0.0636 

1.5117 

7.7469 

-1.4472 

-0.3651 

30.7504 

23.6768 

23.2986 

21.2717 

16.6859 

33.6509 

19.6462 

4 . 9488 

6.5312 

9.8957 

9.9965 

58.2688 

10.7641 

12.8559 

Columns  8 

through  14 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

900.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

400.0000 

10.0000 

10.0000 

10.0000 

30.0000 

30.0000 

30.0000 

30.0000 

80.0000 

200.0000 

500.0000 

30.0000 

60.0000 

90.0000 

300.0000 

14.7739 

12.3204 

-8.4988 

14.6459 

14.4681 

14.2218 

10.5792 

0.0055 

0.0079 

0.0196 

0.0047 

0.0045 

0.0045 

0.0081 

1.0855 

0.6374 

0 . 8890 

1.4480 

1.3821 

1 . 3245 

1 . 1507 

4 . 6071 

6.3695 

7 . 9289 

4.1360 

3.9273 

3.7701 

4 .1144 

-4.1324 

-8.7601 

-12.7254 

-3.5668 

-3.4091 

-3.3334 

-5 . 0978 

0.6303 

3.7443 

5.4575 

0.2898 

0.4899 

0.6914 

2 . 1394 

13.4226 

7.9890 

1 . 1516 

106.1962 

73.7996 

56.8377 

12.2073 

15 . 3371 

27.7931 

31 . 9739 

44.3994 

49.0150 

53.7996 

66.7834 
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Columns  15  through  21 


900.0000 

400.0000 

400.0000 

50.0000 

10.0000 

14.6882 

0.0050 

1.6733 

4.6350 

-5.0731 

0.9066 

185.0244 

61.7791 

Colxomns  i 

900.0000 
400 . 0000 
400 . 0000 
70.0000 
80.0000 
13.2627 
0.0063 
1.8584 
3.4892 
-3.2833 
0.0488 
67.7916 
73.1136 

Columns  : 

900.0000 
400.0000 
360.0000 
10.0000 
7.0000 
14.4258 
0.0051 
1 . 6570 
2.9855 
-0.5888 
-1  7114 
19 . 9502 
22.7010 


900.0000 
400.0000 
400.0000 
50.0000 
40.0000 
14.3058 
0.0048 
1 . 6406 
4.0694 
-4.1116 
0.6769 
111.7432 
65.7842 

through 

900.0000 

400.0000 

400.0000 

70.0000 

200.0000 

11.3093 

0.0094 

1.5995 

2.6174 

-2.1580 

-0.2808 

35.4038 

95.1170 

>  through 

900.0000 
400.0000 
360.0000 
10.0000 
10.0000 
14.4140 
0 . 0051 
1.6234 
3.0660 
-0.7776 
-1.5832 
13 . 6716 
30.6612 


900.0000 
400.0000 
400.0000 
50.0000 
70.0000 
13.9310 
0.0047 
1.6101 
3.5867 
-3.3023 
0.4855 
82.8307 
71 . 8230 

28 

900.0000 
400 . 0000 
400.0000 
70.0000 
500.0000 
7.7383 
0.0248 
1.0749 
4.8419 
-6.5218 
0.6029 
5.8536 
86.7295 

35 

900.0000 
400.0000 
360.0000 
10.0000 
30.0000 
14.2724 
0.0052 
1.4227 
3.5624 
-1.9447 
-0.7848 
27 . 9442 
11.5182 


900.0000 

400.0000 

400.0000 

50.0000 

100.0000 

13.5639 

0.0048 

1.5819 

3.1867 

-2.6455 

0.3325 

66.7305 

78.6520 


900.0000 

400.0000 

400.0000 

90.0000 

30.0000 

13.8442 

0.0064 

2.4179 

3.4468 

-3.0304 

-1.0538 

127.0123 

70.1120 


900.0000 

400.0000 

360.0000 

10.0000 

60.0000 

13.8549 

0.0054 

1.1991 

4.1763 

-3.3961 

0.2304 

19.6615 

14.4526 


900.0000 

400.0000 

400.0000 

50.0000 

400.0000 

10.3139 

0.0136 

1.4235 

3.7394 

-4.4551 

0.9145 

5.0340 

66.6357 


900.0000 

400.0000 

400.0000 

90.0000 

60.0000 

13.1552 

0.0080 

2.2349 

3.7625 

-3.6347 

-0.7201 

91.2764 

75.5835 


900.0000 

400.0000 

360.0000 

10.0000 

90.0000 

13.1911 

0.0056 

1.0682 

4.6334 

-4.4885 

1.0265 

15.0107 

16.9895 


900.0000 

400.0000 

400.0000 

70.0000 

20.0000 

14.3519 

0.0053 

1.9983 

4.1426 

-4.4611 

0.3600 

136.2876 

65.0384 


900.0000 

400.0000 

400.0000 

90.0000 

90.0000 

12.4310 

0.0097 

2.0586 

4.0403 

-4.1715 

-0.4271 

96.0007 

85.1464 


900.0000 
400.0000 
360.0000 
30.0000 
20.0000 
15.1308 
0.0046 
1.5394 
4.0315 
-3.6417 
0.4220 
117 . 4703 
42.4843 


900.0000 

400.0000 

400.0000 

70.0000 

50.0000 

13.7979 

0.0057 

1.9275 

3.7895 

-3.8210 

0.1922 

90.1165 

68.3909 


900.0000 

400.0000 

400.0000 

90.0000 

300.0000 

6.3745 

0.0213 

1.0123 

4.9209 

-6.0402 

0.4859 

23.7164 

91.1735 


900.0000 
400.0000 
360.0000 
30.0000 
50.0000 
14.2001 
0.0043 
1.3707 
3.9955 
-3.8458 
0 . 8102 
79 .2965 
47 . 9708 
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